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ABSTRACT 


1 


A two-dimensional model is developed to describe the charging of thin 
polymer films exposed to a uniform mono-energetic electron beam. The 
study is motivated by observed anomalous behaviour of geosynchronous 
satellites which has been attributed to electrical discharges associated 
with the differential charging of satellite surfaces by magnetospheric 
electrons. 

Electric fields both internal and external to the Irradiated specimen 
are calculated at steady state in order to identify regions of high electri- 
cal stress. Particular emphasis is placed on evaluating the charging 
characteristics near the material's edge. 

The model has been used to identify and quantify the effects of some 
of the experimental parameters notably: beam energy; beam angle of in- 

cidence; beam current density; material thickness; and material width. 
Simulations of the following situations have also been conducted: positive 

or negative precharging over part of the surface; a central gap in the 
material; and a discontinuity in the material's thickness. 
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I . INTRODUCTION 

Gtoaynchronouf iatcllltes in the early aeventieii exhibited ano<oalous 
behaviour which at that time defied explanationt Telemetry indicated 
non-existent status changes » antenna gains changed in the absence of 
commands I tracking and spin stabilization were momentarily lost and 
occasional loss of communications led to total mission failure [Rosen . 1975]. 
Of more than 150 communications and remote sensing satellites now in geo- 
synchronous orbit about 30 have had unexplained and seemingly random 
anomalies. Correlational studies indicated that these failures occurred 
mainly in the midnight-to-dawn local time period and that electron spectro- 
meters in orbit simultaneously registered unusually high fluxes of ener- 
getic electrons. In other words, the problem was linked to an environmental 
effect rather than random component failures [Sharp et al . 1970]. Analysis 
of the ATS-5 satellite plasma detector data by De Forest [1972] suggested 
the possibility that such satellites could be charged by ionospheric electrons 
to potentials greater than 10 kV. Actual confirmation that the charging 
and violent discharging process does occur in geosynchronous orbit was 
established in recent reports from the ATS-5, ATS-6 and SCATHA satellites 
[Stevens . I960; Koons . 1980], 

Synchronous satellite orbits are circular, with a radius of 6,6 earth 
radii. In the midnight-to-dawn sector of such orbits, the satellite en- 
vironment can become filled with a hot plasma during magnetic substorms. 
Electron energies encountered range from 3-200 keV [Hirshberg and Colburn. 
1969] during such intervals and the accumulated satellite charge can there- 
fore result in high spacecraft potentials, The electron current density 
is approximately 0.1 to 10 nA/cm^ [McPherson and Schober . 1976]. The pro- 
ton current density during the substorm is estimated to be only 3% of the 
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fsieccton current dengity [Inouya . 1976], 

* 

the abrupt loss of charge (diacbarge) generatea electromagnetic inter- 
ference which can couple into control systems and temporarily degrade their 
performance. Repetitive discharging has caused some surface damage to such 
materials as optical surfaces, solar cell arrays and thermal control sur- 
faces thus possibly reducing their useful lifetime. The potential gradient 
that develops about the satellite can also perturb normal magnetospheric 
conditions sufficiently to render useless some scientific measurements of 
the environment. 

An understanding of the behaviour of materials under charged particle 
bombardment is essential to the prediction and prevention of the adverse 
effects associated with spacecraft charging. Such an understanding has 
been, sought through ground based tests on the charging and discharging 
characteristics of dielectrics when exposed to uniform monoenergetic elec- 
tron beams [Balmain and Dubois . 1979; Aron and Staskus . 1979; Stevens et al . 
1977; Robinson , 1977; Purvis et al . 1979]. Typical experimental layouts 
employ small (< 100 cm^) planar samples mounted on a conducting substrate. 
Measurements indicate that a large fraction (25 - 50%) of the trapped charge 
can leave the surface in the form of blowoff electrons during the discharge 
event, a phenomenon referred to as "charge cleanoff". Spanning three 
orders of magnitude of exposed-area variation a consistent set of empirical 
scaling laws describes the substrate return current during the discharge; 
the total charge released is proportional to the sample area and the peak 
current is proportional to the sample radius [Balmain and Dubo is. 1979; 
Flanagan et al . 1979]. The appearance of the discharge arcs suggests that 
the discharges start at the dielectric edges or seams and propagate across 
the surface [Stevens . I960]. A second class of discharge, termed 
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"punchthrough", rtltasct a smallar fracCion of charga and if characterlzad 
by cha appaaranca of a pinhola-tlzad punccura in eha matarial. Subaaquent 
dlschargas appaar Co originate in rapid suaeassion from thia aita. 

PaaC modelling of cha laboratory aimulationa haa concancratad on a 
onaodiaanaional daacriptlon of Che charging procafci [Baara at al . 1980; 

Purvia at al . 1977] which raprcaenta the charging procaaa near the centre 
of Cha apaciman. Howavar, cha tendency of the diacharga to ba aaaociatad 
with a wall defined edge auggaata that auch a modal cannot adequataly daacribe 
the conditions of discharge initiation. In addition, data from the SCATHA 
satellite suggests that discharges can occur at differential voltages of 
only 3 kV [Stevens . 1980l. The resulting bulk field in the material cal- 
culated by a one-dimensional theory is not thought to be sufficient to 
cause the dielectric br<$ 4 kdpwn associated with such discharges. Investi- 
gators attempting to understand the discharge process have therefore been 
forced to examine the details of the internal charge distribution and to 
determine multi-dimensional effects which could influence the material’s 
stability [Beers et al . 1981; Stevens . 1981]. 

This work represents an attempt to quantify the higher-dimensional 
effects present in the charging of planar dielectric films exposed to a 
mono-energetic electron beam. Regions of high charge concentration and 
intense electric fields are identified under a variety of charging con- 
ditions in order to isolate factors related to the discharge process. 

The following two chapters (Chapter Z and Chapter 3) are a review of 
the present state of the art in electron beam charging effects. New 
material is introduced in Chapter 4 and the following chapters. 
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2, ELECTRON BEAM CHARGING EFFECTS 
IN DIELECTRICS 

2 . 1 Introductlcig r 

Electron beam effects begin with the interaction between an energetic 
incident electron (1 - 100 keV) and a target atom resulting in a transfer 
of energy. This energy can appear in a number of forms including excita- 
tion of the target atom, vibrational and rotational motion of the atom mani- 
fested as heat, and ionisation of the atom, Some ionizing events transfer 
sufficient energy to create recoil electrons in the material, These can 
initiate further ionizing inters'ictions resulting in an abunwi^nce of free 
electrons which can be moved under an applied field through the usual thermal 
activation process, 

If the primary electron mean free path between collisions is too small 
then its range will not be sufficient to penetrate the material. These 
trapped electrons result in an accumulation of net charge and the development 
of large fields within the dielectric. The remainder of this section is 
devoted to listing some of the material responses which can Influence the 
charge accumulation Ir, the exposed material, 

2.2 Secondary Emission and Backscatter 

Any material exposed to an electron beam is found to emit electrons 
with a broad spectrum of energies, as illustrated in Fig, 2,1 r Harrower . 
1957], The distribution is distinctly bimodal for primary energies great- 
er than 500 eV, Those electrons emitted with lower energy are referred to 
as the secondary emission component while those with higher energy are 
referred to as the backscatter component of the total emission, The rela- 
tive magnitude and characteristic energies depend on the material and the 
^ energy of the primary electrons, 

The high energy of the backscattered electrons is the result of near- 
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tlaitlc collifioni wich th« mattrlal. Th« backscAttar ipactrum hat fcha 
primary anargy at ira naxiaumi Thia high anargy natura iapliat Chat the 
•catcaring procatt could occur over an appraciabla fraction of the pri- 
mariat* range. The backacatter coefficient BS defined at the fraction 
of incident electrona anitcad in the high anargy tail of th$ amitsion 
spectruisi ia uaually found to decraaae tlowly with incraaaing energy. 

For moat polymer dielectrica the following empirical ralationahip holda 


[Wall at al. 1977] 


BS - O.l(KE) 
where BS 


- 0.2 


( 2 . 1 ) 


Jbs^’^i 


backscatter coefficient 


K£ » kinetic energy of primary elec- 
trona (keV) 


The backscatCer coefficient for primaries incident at an angle $ with 
respect to the surface normal of moat metals has been described as 
[Darlington . 1975] 

BS(e) - B(BS(0)/B)‘^°®® (2.2) 

where B p dimensionless material constant 

9 angle of incidence with respect to the 
suirface normal 

if we take B to be unity then this expression also fits Che data for most 
polymers [Wall et al . 1977]. 

The secondary electrons are a direct manifestation of Che ionization 
process in the material. The shape of the secondary emission energy spec- 
trum is largely independent of the primaries' energy. The probability of 
a secondary electron escaping from a given volume of Irradiated material is 
assumed Co decrease exponentially with increasing distance from the material 


secondary backscacter 



500eV BE 


Kinetic energy of einitted electrons 
Fig. 2.1 Generalized emission spectrum. 



KE, KE KE„ 
1 m 2 


Kinetic energy of primary electrons 

Fig. 2.2 Generalized secondary emission yield vs. primary 
energy. 


surface. The decay constant for most polymers varies betwcr^i 0.002 and 
0.005 \m [Budd . 1981]. 

The secondary emission coefficient SE, defined as the number of sec- 
ondaries emitted per incident primary, is a strong function of the primaries 
energy as indicated in Fig. 2.2. The energies where the coefficient is 
equal to unity, K£^ and KE^, are known as the first and second crossover 
points. The mechanisms responsible for the maximum are a minimum ionization 
energy requirement and the Increasing penetration of the primaries at high- 
er energies. Burke [1980] has demonstrated that a universal curve can be 
generated relating the secondary emission coefficient and the primaries' 
energy for most inorganic Insulators. At high energies the secondary 
emission coefficient is given by a simple power law relationship 


SE 


K(KE) 


-0.725 


where KE ■ primaries' energy (keV) 


The constant K is found to decrease as the complexity of the repeating 
unit in the polymer increases. Quantitatively this can be expressed in 
terms of the gram molecular weight M and the number of valence electrons 
in the repeating unit N(l,4,5,6 for U,C,N,0 respectively), Polymers con- 
taining just these elements are governed by the following relationships 
iBurke. 1980]; 


emission coefficient 

K - 

10.64(N/M) - 3.15 

(2.3) 

energy at maximum yield 

KE 

m 

- (K/12.09)°’^® 

(2.4) 

maximum yield 

SE 

m 

- 9.5 KE 

m 

(2.5) 

secondary emission 

SE - 


(2.6) 


where Z * 1.284 KE/KE^ 


POLYMER 


FORMULA 


K 



Teflon 

CF 2 

1.546 

0.303 

Polyethylene 

CH 2 

1.370 

0,283 

Polyvinyalcohol 


1.242 

0.267 

Nylon 

CgHiiON 

1.148 

0.255 

Lucite 


1.115 

0.251 

Polystyrene 

*^8«8 

0.907 

0.223 

Mylar 

‘=10«8°4 

0.847 

0.214 

Kapton 

S2»10°5«2 

0.682 

0.189 


Table 2.1 Emission coefficient K (from Burke [I960]) 
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Fig. 2.3 Experimental results and eq. 2.6 (from Burke [1980]) 
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For KE greater than 1.5 KE„, the following approximation repreients an over- 
estimate of less than 5^ 

SE • K (2.7) 

Values of K for a large number of polymers are presented in Table 2.1. 

In Fig. 2,3 the universal emission curve, represented by eq. 2.6, is pre- 
sented. 

Robinson and Budd [1980] have fitted measurements of total emission 
from Teflon for large primary energies to a single power law term 

BS + SE - (KE2/KE)°'^® (2.8) 

where KE 2 is the second unity crossover energy 
(KE 2 » 2.1 keV for Teflon) 

This figure is slightly larger than that suggested by Burke (eqs. 2.1 and 
2.7) though both predict the same second crossover energy. Budd [1981] 
has also reported measurements of the total emission dependence on the beam 
angle of incidence. For 9 I'ess than 50® it was found to vary as l/cos9 as 
would be consistent with the number of secondaries being proportional to 
the deposition of energy in a fixed-thickness surface layer. Wall et al 
[ 1977 ] suggest an alternative semi-empirical relationship which can be com- 
pared to Robinson's formulation: 

Wall et al SE(9) - SE(0) exp(2(l - cos9)) (2.9) 

Robinson SE(9) - SE(0)/cos9 (2.10) 

For angles in the interval 0® - 78° the functional form suggested by Robinson 
produces a smaller emission yield. 

Surface contamination and surface roughness can also influence secondary 


Second unity emission energy 
KE^ (keV) 


II 
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Normal electric field of surface, E 


Fig. 2.4 Dependence of unity crossover energy on normal 

surface field strength for Teflon (from Robinson 
and Budd [1980]). 



mission [Budd . 1981; Dsvies . 1979]. The possibility of a surface electric 
field dependence also exists due to the relatively low energies of the 
freed secondaries. Robinson and Budd [I980l have made measurements Indi- 
cating that the second unity crossover point is a weak function of the nor- 
mal surface field when accelerating the secondaries away from the surface 
(see Fig. 2. A). Robinson and Nguyto [1979] have also speculated that sur- 
face tangential fields may enhance secondary emission. Balmain [1978] and 
Feder [1976] have used the enhancement of secondary emission due to charge 
storage in a material to image the charge deposited in a dielectric with a 
scanning electron microscope. No quantitative data concerning these effects 
have been found in the literature. 

2.3 Electron Range 

An electron, due to its low mass, is prone to large changes in its 
momentum in any interaction. The statistical nature of these deflections 
results in a very broad deposition profile known as range straggling. The 
number of electrons penetrating to a particular depth in a material does 
not drop sharply at the mean range as is the case of heavier particles. In- 
stead the loss is gradual as shown in Fig. 2.5. 

Various measures of electron penetration have been used including 
mean or average range, extrapolated or practical range and linear range. 

The mean range refers to the distance at which 50 % of the incident electrons 
have been absorbed. The extrapolated range is obtained by projecting the 
slope of the fraction penetrating versus distance curve to a zero 
crossing as indicated in Fig. 2.5. 

The linear range is a calculated quantity based on a formulation of the 
energy loss per unit distance of penetration. The most common approach used 
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in this calculation is to assume a continuously slowing down approximation 
(GSDA range) where the energy loss per unit distance is a constant equal 
to the mean rate of energy loss. The linear range is a measure of the 
length of a possible electron path which can be deflected far from the pro- 
jected Incident trajectory and is therefore usually 25 - 50% larger than 
the extrapolated range. For Teflon the average range is found to be about 
60% of the extrapolated range [Gross et al . 1974]. 

At primary energies less than 50 keV a reasonable representation of the 
extrapolated range is given by the Gledhlll relation [Gross . 1980] : 


log r ■ 
* ex 

where r 

ex 
KE ■ 


-5.1 + 1.358 log KE + 0.21 5 (log KE)2 
- 0.043(log KE)3 

extrapolated range x material density (gm/cm^) 
primary energy (keV) 


( 2 . 11 ) 


From eq. 2.11 for energies less than 10 keV 


avg 


• 0.6 rex 

» 4.8 X lO’^(KE)^’^^® (gm/crn^) 


( 2 . 12 ) 


The electron deposition profile has been calculated using Monte Carlo 
simulations [Beers et al . 1979; Frederickson . 1979; Berkley . 1979] although 
the necessary penetration-energy loss mechanisms for very low energy elec- 
trons are not well established [Beers et al . 1980]. Typical deposition pro- 
files for various primary energies are given in Fig. 2.7 as calculated by 
these codes. 


2.4 Dose Rate 

The dose is defined as the mean energy Imparted by ionizing radiation 
to the material in a volume element divided by the mass in that volume. 





(• IMMIMI -•/• 0I> - ■»/• -01) • 



Fig. 2.7 CHARGF DEPOSITION PROFILE IN TEFLON 
FOR NORMALLY INCIDENT ELECTRONS 

(from Beers etal [1977]) . 



We shall adopc as a measure of dose the "rad" (1 rad ■ 10“2 Joules/kg) , The 
average absorbed dose rate for non-penetrating, normally incident electron 
beam is [Gross . 1980]; 

dose rate - 6 • I^KE k 10®/(a“ 6 ) (rads/sac) (2,13) 

where D > dose (rads) 

Ij. ■ Incident current measured to the material 
surface a (A) 

KE • primary electron energy (keV) 
a ■ surface area (cm^) 
r ■ average range (cm) 

S ■ material density (gm/cm^) 

Assuming a constant current to the surface eq. 2.13 can be modified for non- 
normal incidence by dividing by the new irradiated volume aTcosS, Instead 
of ar: 

D(e) - D(0)/cOse (2.14) 

If the range-energy (eq. 2.12) and the dose-energy (eq. 2.13) relations are 
combined we find that the mean dose rate decreases as the energy of the 
primaries increases. 

Details of the actual deposition profile can again only be calculated 
using Monte Carlo techniques [Beers et al . 1979; Frederickson . 1979]. Typi- 
cal curves show a broad maximum near ths mean range due to beam straggling. 

2.5 Conductivity 

The conductivity g of most polymers under low electric fields, in the 
dark and at room temperature is very small making it often difficult to 
establish accurately. For Teflon FEP the estimates vary considerably: 
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< 10*^^ Cncrn)”^ [S*s«1t . 1974] 

3.3 X lO"^® (ncm)”^ [Beers et al . 1979] 

(0cm)"^ [Kinney . 1957] 

10“^^ - lO”^^ (ncm)“^ [Frederlckson . 19791 

For Kepton H there has been less work published although a figure between 
10**^^ and 1Q~^^ (ncm)**^ has been established from s series of measurements 
made by Adorno and Nanevicz [1977]. In general it has been observed that 
Teflon is a better Insulator than Kapton by &bout two orders of magnitude. 

2.5.1 Conductivity: Temperature Dependence 

The conductivity of most insulators displays a temperature 
dependence that can be described by the following empirical expression [Wall 
et al . 1977]: 

- A exp(-W/kT) (2.15) 

where T ■ absolute temperature >, 300°K 
A,W ■ material constants 

In the limit of zero applied fields Hanscombe and Calderwood [1973] found 
W to be 1.55 eV for Kapton. Amborskl [1963] found W to be about 1.0 eV for 
Kapton under a field strength of 5 x 10** V/cm. Using data from Adamo and 
Nanevicz [1977] we find W - 0.77 eV for Teflon FEP under fields of 10^ V/cm. 

Under normal spacecraft operating conditions temperature ranges can be 
expected from 0°C to 120°C [Rosato . 1968; S treed and Arvensen . 1967]. Using 
W equal to 1.0 eV we find that over this temperature range the magnitude of 
the dark conductivity could Increase by a factor of 5 x 10®. 

2.5.2 Conductivity: Photoconductivity 

The effect of light on the conductivity of a variety of materials 
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ha» b««n r«port«d brlafly bv Adano and Nanevlcz (1977). The moit eeneitlve 
material listed is Kapton which when illuminated by a Xenon lamp at an in- 
tensity of one sun displayed an increase of up to 4 orders of magnitude in 
conductivity. There was no attempt to isolate temperature effects or spec- 
tral response. Ho further work has been found in the literature. 

2.5.3 Conductivity: Field Dependence 

The conductivity of most polymers is a non-linear function of 

applied electric field. In particular there exists a range of high fields 

1/2 

where it increases as exp(E ). A functional form describing this be- 
haviour is [Adamec and Calderwood. 1975] 

1/2 

g/go - 121 ) (.^sinh(bE/T)) (2.16) 

where a,b ■ material constants 

T ■ absolute temperature (°K) 

E ■ applied electric field 

For fields less than 10® V/cm the second term in eq. 2.16 is near unity. 

The value of a/T at 22°C was found to be 4.68 X 10“^ (cm/V)^^^ and 

-3 

4.03 X 10 (cm/V) for Kapton and Teflon respectively as indicated in Figs. 
2.9 and 2.10. 

2.5.4 Conductivity; Radiation-Induced Conductivity 

The charge carriers generated by the ionizing radiation give 
rise to an enhanced conductivity. A power law relationship between the 
equilibrium radiation-induced conductivity and the dose rate has the follow- 
ing form 

*rlc " 

where A.G ■ material constants 
0.5 <. 4 <. 1.0 

D ■ dose rate 


Conductivity, g(ncm) Conductivity, g(fic«^ 


y>] 


5x10 - (2+cofh(*,68xli0 ' E^))/3 


• t 

t •• 


, * *• ttoparaturc 

: ‘ 22<5c 

• data pointg from 

Adano and Nanavicz (1977] 


Eleetxic field, E (V/cm) 

Fig, 2.8 Field induced conductivity in Kapton. 



Electric field, E (V/cm) 

Fi*. 2.9 Field induced conductivity in Teflon. 
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The value of h la determined by the energy diatribution of trapa in the 
material; a uniform diatribution givea ^ equal to l.O, while a aingle ttap 
level givea L equal to 0.5 [Qroaa . 1980]. It la often aaaumed that the 
radiation induced conductivity la proporational to the done rate ( A ■ 1) 
tBeera . 1979; Mall et al . 1977]. Table 2.2 Hats the range of values 
found in the literature for G baaed on this aasumption. 

Groat et al [1974] report that meaaurementa made at low dose rates 
tend to yield lower values of t than those made at high rates. Data in 
Figi 2.11 produce a power law fit for Teflon FEP at an average dose rate 
of 2 X 10** rad/sec; A « 0.7, G - 1.7 x lO"^^ (Gcm)"^. In fig. 2.12 this 
is plotted together with the upper and lower limits suggested by Table 2.2. 
The choice of A ■ 0.7 seems to fit the bounds quite well over the dose 
rate interval of Interest. Fitting a similar curve for Kapton we find 
A o 0.7, G - 1.7 X 10"^^ (ncm)”^ 

A delayed component of radiation induced conductivity persists after 
the radiation has been removed for a period of seconds due to the finite 
lifetime of the free charge carriers. A functional form suggested by Gross 
[1980] is; 

8ric “ 8’/(l+bt) (2.18) 

where g* ■ 1/3 of the radiation induced 

conductivity at the end of irradiation 

t ■ time after end of irradiation (sec) 

b ■ constant approximately unity (sec) ^ 
[Beers et al . 1980] 

The neu conductivity of the sample is calculated as the sum of the 
radiation induced and dark conductivities, 
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Fig. 2.11 Functional form of radiation induced conductivity 
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Table 2.2 Upper and lower limits of G. Assumed functional 
form g - g 6 (from Wall et al [1977]) . 
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2.6 Optical Effects of Irradiation 

Irradiation by an electron beam of sufficient energy can cause any of 
the following optical effects [van Lint et al . 1980]; 

luminesisence 
Cerenkov radiation 

introduction and/or bleaching of optical bands 
changes in optical scattering 
changes in ttie index of refraction 

Luminescence for the energy range of interest is easily detected for most 
of the materials under study here. The other effects are small except at 
very large doses or high energies. 

2.7 Radiation Chemistry 

The subject of radiation chemistry is too broad to be treated in detail 
at this level. For a more detailed discussion the reader is referred to the 
literature [ Ausloos. 1968], In polymers ionizing radiation produces free 
radicals that can effect chemical changes falling under two classifications; 
cross-linking and scission. Cross-linking refers to the formation of bonds 
between adjacent polymer chains. Scission refers to the breaking of a chain. 

Frederickson [1979] claims that Teflon starts to degrade when the 
accumulated dose exceeds 10^ rads, while other dielectrics require more than 
10® rads. Wall et al [1977] suggest that variations in the reported values 
for the radiation induced conductivity constant G may be caused by radi- 
ation damage. Gross et al [1974a] report that large doses of the order of 
50 X 10® rads can cause a reduction in G by at least an order of magnitude* 

A 10® rad dose can be delivered to a 0,1 ym deep layer by a 2 keV electron 
beam when the current density - time product is 65 nC/cm^. A 10 nA/cm^ beam 
can therefore deliver such a dose in 6.5 sec. 
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2t8 Electrical Breakdown and Surface Discharge 

Electron Irradiation of Insulators can produce an accumulation of charge 
sufficient to cause dielectric breakdown which Is accompanied by the removal 
of a large fraction of the stored charge. The first published description 
by Gross [1957] Indicated that the discharge was localized at a depth 
roughly equal to the mean range of the primary electrons. Further experi- 
mental work has evolved to help quantify the discharge process: Balmain and 

Hlrt . 1980; Balmain and Dubois . 1979; Balmain . 1978; Purvis et al . 1977; 

Roche and Purvis . 1979; Yadlowsky et al . 1980; Aron and Staskus . 1979; 
Stevens . 1980; Bosma and Levadou . 1979. The actual discharge mechanism, 
however, remains largely unexplained. There is evidence that defects in the 
material may play a large role In the breakdown process [ Davies . 1977]. 

The breakdown threshold of most dielectrics has been found In the range 
1 X 10® - 5 X 10® V/cm [Beers et al . 1979; Davies . 1977; Frederlckson . 1979]. 
Reports also suggest a close association of the discharge with the edge of 
the tjxposed material [ Beers et al . 1980; Stevens et al . 1978; Balmain . 1978; 
Robinson . 1977]. Stevens [1980] concludes that the data collected from the 
SCATHA satellite indicate breakdowns in orbit are occurring despite the com- 
paratively low (2 - 4 keV) differential voltages that have been observed, 
implying that the surface potential may not be the deciding factor In the 
initiation of the discharge event. In turn prompting questions about the In- 
ternal charge distribution and multl-dimenslonal effects. 
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3. ONE-DIMENSIONAL MODELLING 

A onc-disien«lonal study of ths charge build up on planar dielectrics 
exposed to a uniform monoenergetlc electron beam has the advantage of sim- 
plicity over a multi-dimensional treatment. Such models can be used to 
check parameter sensitivity, to determine main factors governing the 
charging process, and to serve as a standard to which the results of more 
complex analyses can be compared. 

The analysis can be separated into an external and an internal charging 
model. The external model is concerned with the net accumulation of charge 
in the material and deals with the measurable quantities of surface po- 
tential and total charge. The Internal model is concerned with the distri- 
bution of charge within the material and the resultant internal electric 
fields. The work in this chapter follows that outlined by Beers et al [1979]. 

3.1 External Model Development 

The net accumulation of charge on a planar dielectric exposed to a uni- 
form mono-energetic electron beam closely resembles the charging of a simple 
resistor-capacitor network. The resistive elements represents the mechanisms 
of charge loss from the surface and the capacitive element corresponds to 
the physical geometry shown in Fig. 3.1. To determine the capacitance per 
unit area we assume that the structures are infinite in a x-z plane thus 
obtaining a one-dimensional differential equation derived from Gauss' Law 
describing the relationship between the electric field and the volume charge 
density 

3Ey/3y ■ -q(y)/e (3.1) 

where Ey is the electric field normal to the 
surface 

q is the electron volume charge density 
(C/m^, positive for electrons) 
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note; since the electron is the principal 
charge carrier the currents and charge den- 
sities will be referenced to the movement 
or accumulation of negative charge, to conform 
with the major part of the relevant literature. 

We shall assume that the space charge density due to the beam current is 
negligible. Thus in the interval (-L,0") the electric field is constant. 
The grounding of the plates at y ■ -L and y ■ d imposes the following con- 
straints on the electric field. 


I'd 


E(y)dy - 0 


^-L 


or 


E(y)dy - -EqL 


(3.2) 


0 


where <■ magnitude of the electric field for 

-L < y < O” 


From eq. 3.1 we find 

E(s) - ECO'*') 


q(y)/eo^rdy » s > 0+ 


^d 


E(s)ds - E(0+)d 


ds 


q(y)/eoEr‘Jy 


(3.3) 


-EqL - E(0''')d 


(S 


ds 


q(y)/eoMy 


Using the condition of continuity for the normal component of the displace- 
ment density we find E(0''’) - Eg/ey. Integrating the right hand side of eq. 
3.3 by parts we find 


-V 


(1 + d/(LE^))eoer 


(3.4) 
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whtre V ■ 


P ■ 


y " 


*d 

E(y)dy the potential at the dielectric 
''0 surface 

'd 

q(y)dy effective charge per unit area 
'd 

yq(y)dy/p the charge centroid 


For thick samples where d » y the external variables (voltage and 
charge density) are largely Independent of the internal charge distribution. 
The range-energy relationship illustrated in Fig^ 2.6 suggests that unless 
considerable low energy excess charge transport occurs after deposition 
the approximation that d » y will be valid for samples thicker than 25 pm 
and beam energies less than 20 keV. If we also assume that the surface 
capacitance to the ground at y ■ -L is negligible compared to the capaci- 
tance to the ground at y - d then eq. 3. A reduces to 


p - -V G (3.5) 

where C ■ eoCr/d the capacitance per 

u^iit area 


Restricting the possible charge transfer mechanisms to the incident 
beam, conduction current, secondary emission and backscatter we find: 


r 



where 

net i bs se c 

Ji ■ Incident beam current density 
J|,g* backscatter current density 

secondary emission current density 
Jg p conduction current density 


(3.6) 


(3.7) 
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Combining cqs. 3.6 and 3.5 

fV 

dV/J_.tW - -(t - to)/C (3,8) 

jy nc 

o 

Equation 3.8 corresponds to that previously used by both Purvis et al [1977] 
Beers et al [1979] to describe the net charge accumulation in the 
dielectric. 

It is implied in eq. 3.8 that the net current density is a function of 
the surface voltage. This voltage dependence is due in part to the secondary 
emission and backscatter currents which tend to increase with decreasing 
primary kinetic energy. The decrease in the impact energy is directly related 
to the potential barrier, i.e. surface potential, that must be overcome to 
reach the surface. By performing an energy balance we find: 

KE + |eVl - BE (3*9) 

where KE - electron kinetic energy at impact 
V ■ surface potential 
e - charge on an electron 
BE ■ the total beam energy (a constant) 

In general the net current density is a non-linear function of the surface 
potential and numerical integration must be used to solve eq. 3.8. 

The model does not incorporate the following effects: 

- surface conductivity 

- radiation hardening of the secondary emission coefficient 

- field dependence of the secondary emission coefficient 

- temperature effects 

- photoconductivity 

- photoemissipn 

- material degradation 

Any of the above could be incorporated once a reasonable data base had been 
established. The development of eq, 3.8 involved several assumptions which 
are inherent in this discussion. These include; 


- no multidinonsional effects 

- no spacecharge effects due to current densities 

- quasi-static potential (potential does not apprsciably 
change during electron transit tine) 

- thick-film, short range approximation 

Two materials, Teflon and Kapton, will be studied in detail. Teflon has 
a relative permltivity of approximately 2.25 and a density of 2.2 gm/cm^ 
permittivity whereas Kapton has a relative permltivity of approximately 3.6 
and a density permittivity of 1.42 gm/cm^ (Beers_et_al, 1980]. The expressions 
for the backscatter and secondary emission coefficients will be chosen to agree 
with those proposed by Burke [1980] and Wall et al [1977] (i.e. eqs. 2.1, 

2.2, 2.7 and 2.9). 

3.2 External Model Results 

3.2.1 Steady State Surface Potential 

A steady state is achieved by reducing the net current to the 
surface to zero, suggesting from eq. 3.7 that 

*^i " “^bs *^se 
or equivalently 

J^d - BS - SE) - (3.10) 

For a normally incident beam, eqs» 2.1 and 2.7 can be substituted for BS and 
SE using eq. 3.9 to express the kinetic energy at the surface in terms of the 
surface potential and the beam energy. Assuming the steady state kinetic 
energy is greater than 1.5 KEj^ (a necessary condition for eq. 2.7 to be 
valid) we obtain: 

J^d - O.KBE - leVl)"°*^ - K(BE - | evj )"°*'^^^) 

(3.11) 

- lv|/d • g(v/d) I 

where g ■ material-dependent conductivity 
K ■ material constant 
0.68 keV for Kapton 
1.55 keV for Teflon 
V ■ surface potential (kV) 


30 


The bulk conductivity of Teflon is assumed to be a constant go in order to 
simplify the following analysis. We shall relax this constraint allowing a 
field-dependent conductivity for the analysis of Kapton. 

The solution of eq. 3.11 for Teflon is presented in Fig. 3.2. Two 
distinct regions exist depending on the ratio of conduction current to beam 
current. If the beam current dominates then the steady state surface poten- 
tial will be proportional to the difference between the beam energy and the 
second unity crossover energy for total emission (BS -f SE) . If the currents 
are comparable, as is the case for low beam currents, the surface potential 
is proportional to the beam current. The asymptotic behaviour is given by: 

if J^/Jc *1, jv| - Jid/ggCl - SE - BS) evaluated for KE - BE (3.12) 

if Ji/Jc » 1, |V| - (BE - KE2)/le| (3.13) 

where KE 2 ■ second unity crossover energy for 
total emission 
0,69 keV for Kapton 
2.10 keV for Teflon 


The critical value of J^d/g^ separating the conduction-limited and emission- 
limited regions is obtained by equating eqs. 3.12 and 3.13. For Teflon we 
f ind : 


So) critical 


(BE - 2.1)/|e 


(1 - 1.55BE-0*^25 _ o.IBE-0-2) 

where BE ■ beam energy (keV) 


(3.U) 


Using the total emission coefficient of eq. 2.8 suggested by Robinson and Budd 
[1980] results in a slight downward shift in the conduct IvityTlimited region 
of the curves. 

The introuduction of a field dependent conductivity as suggested by eq. 


Fig* 3.2 EqulllbcluM surface potential, as a function of the charging paraa«terS J 
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2.16 foe KapCon dots not afftet tha asymptotic analysis pravioualy discussad. 
The transition batwaan thasa ragions howavar, is much expanded as indicated 
in Fig. 3.2 for Kapton. Tha affect is most significant for high beam 
energies and thin samples. 

Interpreting Fig. 3.2, we find that for typical spacecraft charging 

-18 —1 

parameters (material thickness 25 ym, bulk conductivity 10 (Hem) ) the 
critical current density for Teflon is 0.013 nA/cm^. For densities greater 
than this the equilibrium surface voltage and charge will be practically in- 
dependent of the beam current. Balmain and Hirt [1980] have found that the 
total charge released from Teflon during a discharge event displays this 
same independence for typical charging conditions (current density O.A - 100 
nA/cm^, beam energy 20 keV, material thickness 30 ym) as would be expected 
if the specimen is emission-limited, Balmain and Dubois [1979] have also 
found that the total charge released is proportional to the area irradiated 
(see Fig. 3.3). Figure 3.3 indicates that approximately 25% of the net 
available charge at equilibrium in Teflon is released per discharge. For 
Kapton with a beam current density of 80 nA/cm , the charge released is 
reduced to 10%, 


In Fig, 3,4 the equilibrium surface potential for Kapton has been graphed 
in the transition region (assuming gg • 5 x 10 as suggested in Fig. 2.9 
and a beam energy - 20 keV) . The span of current densities, 0.1 - 100.0 nA/cm^, 
covers the typical range employed in laboratory simulations. A discharge 
dependence on the incident current density is therefore expected and Indeed 

Balmain and Hirt [1980] have reported a total released charge dependence of 

n 21 

di as Indicated in Fig. 3.5. The total available charge, under similar 
charging conditions as indicated in Fig. 3.4 has a current density dependence 
of 

3,2,2 Charging Dynamics 

The development of the stored charge is governed by the current 


Ckmttfi Q (irf:) 
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Fig. 3.3 Th« variation of ralaaaatf eharga Q • | X dt with dlalaetrlc apaelaan araa. 
(from Balnaln and Dubois [1979]) . 


net charge available at equilibrium 





Fig. 3.4 Equlllbrivun surface potential and charge density 
vs. current density for Kapton. 
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balance to the surface (eq. 3 #6) which is dependent upon the surface voltage 
through the effective capacitance (eq. 3.5), these tv;o relations giving: 


Jnet(V) - -C dV/dt 


(3.15) 


For Teflon the incident current is assumed to be sufficiently large so that 
the charging process is emission-limited. Equation 3 >15 can then be written 


as follows: 


tJ^d/eoCr 


dV 

(1 - BS - SE) 


(3.16) 


Numerical integration of eq. 3.16 for the energy-emission relationships 
suggested by Burke [1980) (eq. 2.1 and 2.7) and Robinson and Budd [198Q] 

(eq. 2.8) produces the curves in Fig. 3.6. Both relationships yield reason- 
ably consistent results. The experimental results of Purvis et al [1977] 
are in general much faster, but show the same qualitative behaviour. A more 
detailed comparison will be made in Chapters 4 and 6. Good agreement is 
obtained with the model developed by Beers et al [i960]. 

The voltage and charging dynamics as suggested in Fig. 3.6 can be matched 
to an arbitrary beam current density or material thickness by a simple scaling 
of the time axis. 


3.3 Internal Model Development 


The development of intense internal fields was first proposed by 
Meulenburg [1976]. It was suggested that two oppositely charged layers could 
develop in the irradiated material. As discussed in Chapter 2 secondary 
emission is a process that is essentially a surface phenomenon whereas the 
electron beam can pass through the surface with little attenuation resulting 
in a depletion of negative charge at the surface as Illustrated in Fig. 3.7. 

A second negative charge layer was proposed at the mean electron range. If 






Fig. 3.7 Internal charging toodel geometry 
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the conductivity between these two layers, being dominated by the radiation 
induced component, is too small then breakdown fields could develop. Such 
a charge differential could continue to grow long after an external equilib- 
rium has been established. These ideas have been refined in later work to 
include the deta:|.ls of electron deposition and beam induced conductivity 
profiles [ Frederlckson. 1979; Beers et al. 1979, Berkley. 1979]. 

The Internal charging problem can be formulated from Gauss' law and the 
equation of continuity: 

9E/3y - -q/c (3.17) 

3q/3t - 3Jt;/3y (3.18) 

where q ■ electron volume charge density 
(C/m3) 

total electron current density in 
the material in the neg. y direction 

Following Beers et al [1979] we can isolate the electric field by combining 
eqs. 3.17 and 3.18. 

32E/3y3t - -1 /e 3J^/3y 

integrating with 

respect to y 3E/9t - -l/e(Jt(y,t) - Jo(t)) (3.19) 

where - Jp(y,t) + g(y,t) E(y,t) 

J ■ primary beam current as determined by 

the electron deposition profile including backscatter 

g * conductivity as determined by the local 
dose rate 

J_(t)« net external current density 
« J^(l-BS-SE) 

Both the primary beam current profile and the conductivity profile are 
influenced by the kinetic energy of the incident electrons which in turn is 
fixed by the surface potential and the beam energy (eq. 3.9). Equation 3.19 
is therefore coupled to the surface potential described by eq. 3.16 

The effect of covering the front surf ace with a thin (relative to the 
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range of the electrons in the covering) grounded metallized layer has been 
studied extensively [Gross et al . 1980; Oliveira and Gross . 1975; Gross et al . 
1974]. A similar analysis is applicable to spacecraft charging conditions 
when photoemission from the illuminated specimen prevents the material from 
charging [Beers , 1979]. The steady state field under such conditions can be 
easily determined by requiring 3q/3t - 0. Using eq. 3.18 we find the equiva- 
lent condition 3J^/3y ■ 0 suggesting that 

Jt - g(y)E(y) + Jp(y) - J (3.20) 

where J is some constant to be determined 


The condition that both faces of the dielectric are at ground potential is 
used to find J which then is used to detemine the internal bulk field v/here 

rd 


J = 0. 
P 


E dy 


J ■ 


Jo V® ■'y 


(3.21) 


^0 


1/g dy 


3.4 Internal Model Results 
3,4,1 Controlled Beam 

By adjusting the beam energy, the kinetic energy at the surface 
can be fixed making Jp and Jg in eq. 3.19 constants in time. The resulting 
equation is easily solved: 


E(y,t) - “(Jn(y) - JJ/gCy) 


(3.22) 


+ (E(y,0) + (Jp - Jg)/g(y)) exp(-g(y)t/e) 

The final solution is independent of the initial field in the material at the 


trii «lifi>i iMH ■)! I II ifimn I f r ti I 
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start of controlled irradiation. 

Most spacecraft charging experiments are not performed with beam control) 

however eq. 3.22 can be applied to a steady state analysis for which Jq Is 

constant and approximately equal to zero. For Teflon this occurs at an 

incident kinetic energy of 2.1 keV which, using eqs. 2.13, 2.11 and 2.17, 

-13 -1 

implies an average radiation induced conductivity of 1.5 X 10 (^cm) in 

the front Irradiated volume for a beam density of 1 nA/cm^. The resulting 
time constant is approximately 1.4 sec suggesting that those fields quickly 
assume a steady state value independent of the earlier charging history: 

E(y) - -Jp(y)/g(y) y < x <3.23) 

where r ■ average electron range 

In the non-lrradiated bulk region the dark conductivity is approximately 

—18 —1 5 

10 (Qcm) giving a time constant of 2.2 x 10 sec (s 61 hrs) . The bulk 

field is essentially frozen at its initial value 

E(y) - E(0,y) y » r 

Under emission-limited charging conditions we have found 

E(y) » (BE - KE 2 )/(d - 7 ) y » 7 (3.24) 

For a non-charging beam (Jq = 0) Beers et al have calculated the steady state 
field profile based on eq. 3.23 and illustrated in Fig. 3.8. 

3.4.2 General Front Face Field 

A unique front face field, first suggested by Beers et al [1979] 
for any beam energy, current density and material thickness, is a consequence 
of the assumptions inherent in the model. The magnitude of the internal 
field at the front face can be calculated at equilibrium as follows: 


Electric Field CIO V/m) 



Fig. 3.8 SATURATION ELECTRIC FIELD IN TEFLON 

FOR NON-CHARGING NORMALLY INCIDENT BEAM 

(from Beers et al,[1979j) 


MATERIAL 

EQUILIBRIUM FRONT FACE FIELD 
(V/cm) 

Teflon 

0.11 - 5.5x10^ 

Mylar 

52.0 - ei.xio^ 

Kapton 

1,8 - 92.xl0^ 


. 

Table 3.1 Range in equilibrium front face fields for the fj 

values of radiation induced conductivity coefficients j 

in Table 2.2 (from Beers et al, [1979]) 
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from eq# 3.19 

ac Che surface 
and In general 

thus 

at equilibrium 


3E/3C 


-l/c(Jp -»■ gE - Jp) 


Jp - Ji(l - BS) 

Jq - - BS - SE) 

e 3E££/3C + gE^£ ■ *• J^SE 
Efi - -j£ SE/g 
where E££ * front face field 


Using eqs. 2.7, 2.12, 2.13, 2.17 and 3.26 we find 


j£ 


1-4 


K KE, 


(0.3584-0.725) 


i‘='ffi rrr 

G(2.08 X 10^3)^ 

where for Teflon KE2 ■ 2.1 keV 


G 

4 

K 


1.7 X 10 
0.7 
1.58 


-16 . 


-1 


(3.25) 


(3.26) 


(3.27) 


If the assumption is made that 4 ■ 1 it is found that E££ is fixed only 
by material constants. Table 3.1 lists the computed bounds for Che front 
face field using detailed dose profile calculations and the range of radi- 
ation-induced conductivity coefficients in Table 2.2 (assuming 4*1) [Beers 
et al . 1979]. Using an average uniform dose based on Che equations presented 
in Chapter 2 the magnitude is bounded between 5.7 x 10 - 2.8 X 10 V/cm for 

Teflon. Using eq. 3.27 with 4 ■ 0.7 and a current density of 1 nA/cm^ the 

3 

equilibrium field is 6.1 x 10 V/cm with a current density dependence of 
(Ji/l nA/cm^)°*^. 

It should be noted that all of these estimates are less than the expected 
breakdown threshold, the largest occurring in Kapton. The discrepancy in the 
results derived here and those obtained by Beers et al [1979] appear to be 
due to the assumed range-energy relationship (eq. 2.11) at the low energy end. 


3 #4.3 Grounded FronC Face Reeulta 

The eceady state field can be determined from eq. 3.20 once the 
primary electron current has been established. For thick samples we find 
from eq, 3.21 that 

Jp(y) » 7 for y < 7 

thus E(y) ■ -Jp(y)/g(y) (3.28) 

where g(y) is dominated by the radiation 
Induced conductivity 

The front face field can be calculated as in section 3.5.2. Using eq. 2.1 
for the backscatter coefficient we find at the front surface: 

Jp - J^(l - BS) • 0.9 Ji 

therefore at steady state l^ffl * '^ 1^8 

- (4.8 K 10"^^)^KE°‘^^®^J^^^-^^0.9/G (3.29) 

Whereas in eq. 3.27 the value of KE was fixed at the second unity emission 
energy, the value of KE in eq. 3.29 is equal to the beam energy since the 
front surface is held at ground potential. 

For ■ 1 nA/cm^. KE « 20 keV and 4 ■ 0.7 the front face field is 
equal to 1.1 x 10 V/cm (roughly twice that found for a floating front face). 
Mote there is a weak dependence on both the beam current density (Jj,/1 nA/cm^) 
and the beam energy (BE/20 keV)^'^^. 

In the non- irradiated bulk of the material the average field strength 
will be reduced by a factor r/d from that found in the irradiated volume. 

For a 50 ym thick sample under the above conditions the bulk field should be 

3 161 
approximately 1,7 x 10'^ V/cm with a beam energy dependence of (BE/20 keV) * . 
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To lunmarize the grounded front turf ace result t^ the front face field 
It increased by roughly a factor of 2 while the bulk field la decreased by 
a factor of 2 x 10^ when compared to the sample calculations made for a 
floating front surface. Beers et al [1979] have calculated an upper limit 
for the front face fields for A ■ 1 as seen in Fig, 3,9, The results are 
independent of current density and use a detailed dose profile calculation, 
Kute for both Mylar and Kapton the upper estimates are comparable to the 
material breakdown threshold at beam energies greater than 5 keV, 

3.4.4 Experimental and Numerical Results 

Experiments performed by Gross et al [1977] have addressed the 
question of internal charge deposition and migration. The charge centroid 
was found to be equal to the average electron range at the start of irradi- 
ation and it subsequently approached the extrapolated range if time was 
allowed for internal charge migration. As the net stored charge increased 
a corresponding increase in the final charge centroid depth was noted, re- 
flecting the higher internal field. This dependence is weak being approxi- 
mately proportional to the logarithm of the bulk field. 

The charge centroid was found to be sensitive to beam current densities 
greater than 5 nA/cm^ decreasing with increasing current density while the 
net injected charge was held constant. Such behaviour supports the use of 
a non-linear radiation induced conductivity as outlined in section 2.5.4. 

Beers et al [1980] have developed the necessary computer codes to solve 
munerically the complete internal charging history for a one dimensional 
model (eq. 3.19) and have presented steady state results for both Kapton and 
Teflon (Figs. 3.10 and 3.11 respectively). The extremes in the internal 
electric field appear at the front face of the material or in its bulk. The 
behaviour of both of these quantities have been discussed in the preceding 
sections. 
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Fig. 3.9 Front face field at equilibrium for both surfaces held 
at ground potential. Minimum value of G from Table 2.2 
(from Beers et al . [19791) 
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Fig, 3,10 Equilibrium electric field profile in Kapton. 
(from Beers et al, [1980]) 
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Fig. 3.11 Equilibrium electric field profile in Teflon, 
(from Beers et al, [1979]) 
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4. TWO-DIMENSIONAL SIMUUTION DESCRIPTION 

A fully chrct-diacniional dcicriptlon such as the NASA charging ana- 
lyser program NASCAP [Katz at al . 1979] can dascrlba the charging of an 
arbitrary body exposed to a general elcctron/ion environment. Unfortunate- 
ly a penalty must be paid for such flexibility in terms of cost and com- 
plexity . 

In this study our attention has focused on structures that can be com- 
posed of long, Uniform parallel dielectric strips above an infinite con- 
ducting ground plane. Under such conditions the accumulation of charge can 
be described in only two dimensions thus achieving considerable savings in 
both cost and effort while isolating the effect of the metal-dielectric edge 
in Che charging process. The two-dimensional model described in this section 
is a direct extension of the one-dimensional work described in Chapter 3 to 
include those effects associated with the deflection of the Incident beam by 
the developing charge distribution. 

The basic geometry, symbol definition and coordinate system are given 
in Fig. 4.1. The results obtained, although idealized, should still be 
applicable to the central portion of any long uniform strip. Many labora- 
tory simulations of spacecraft charging [Balmain and Hirt , 1980; Purvis et al . 
1979; Yadlowskv et al . 198Q] have been performed that can be adequately re- 
presented in this manner. 

4.1 External Charging Model 

The external charging model specifying the accumulation of net charge 
in the dielectric as a function of time is described by the equation of 


charge continuity 

8q/9t ■ 

* "^net 

(4.1) 


where q ■ 

volume charge density 





As was indicated In Section 3.1 the voltage and external eldctrlc fields 
are relatively insensitive to the internal charge distribution for thick 
samples. For the purposes of the external model we shall consider the 
charge to be deposited at a constant height d, equal to the material thick- 
ness. above the ground plane. Using eq. 3.7 for the net current density, we 
can write eq. 4.1 in the following difference form: 

p(z, T+t) - (Ji - Jbs - Jse - 

where all currents refer to the component in the 
negative y direction (the question of tan- 
gential conduction current effects is 
addressed in Section 4.3) 

p ■ net surface charge density 

In general all of the currents are functions of time and position z because 
they are determined by the growing charge distribution. The problem is 
therefore to relate the net current at the surface to the net charge residing 
on the surface. 

A solution outline is gi\'en in Fig. 4.2 in the form of a flowchart. 
Details on each stage of the chart are given in subsequent subsections. 

4.1.1 Derivation of the Electric Field from an Arbitrary Surface Charge 
Distribution 

The charge layer responsible for the unknown field is assumed to 
be located in the dielectric at a constant height above an infinite ground 
plane. In general, the effect of the dielectric is to reduce the effective 
charge layer due to the, polarization of the material. From Jordan and Balmain 
[1968] if the material is homogeneous the potential at a given observation 
point 7 (see Fig. 4.3) due to this polarization is given by 


50 


Fig. 4.2 Flow Chart: External Charging Model 


- set Initial charge distribution to zero 
>■ input beam characteristics 

- input material characteristics 


> given the charge distribution calculate 
the electric field (sec. 4.1.1) 


- given the electric field calculate how 
the defined beam is deflected (sec, 4.1.2) 


y 

- given the beam deflection at the charge 
layer, calculate the local impact para- 
meters; 

kinetic energy 
beam angle 

beam density (sec. 4.1.2) 




- given the Impact parameters and electric 
field calculate the net deposition profile 
across the material (sec. 4.1.3) 


- assuming the given net deposition profile 
is static for the clock period increment, 
calculate the charge distribution profile 
as described by eq. 4.2 


- check if deposition profile is negligible 
relative to accumulated charge 


if no * 


if yes, output final charge 
distribution 
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where pp ■ 'p • n 

P ■ polarization 


Equation 4.3 suggests that the dielectric can be replaced by a surface charge 
equal to Pp. In addition we can replace the ground plane by the appropriate 
image charges to obtain a simple volume charge distribution in free space. 

The final potential is calculated using the resulting charge distribution in 
free space and the familiar Green's function: 


itotal 


p/4ttEoR dV 


(4.4) 


where p ■ equivalent charge in free space 


This sequence of steps is illustrated in Fig. 4.4. Comparing Figs. 
4.4(a) and 4.4(d) we find 


So(Ei - - -(Pp + Pg) 

(4.5) 

“ ^r®2) ■ “Pg 

therefore 

(1 - E,/EJ 

(Pp + Ps)/ps - Tg ; - j-^e ~ T 

where pp Hh Pg •• equivalent charge in free space 

pg « actual net free charge in the dielectric 


Unfortunatelyj Ej^ and E^ cannot be determined until the equivalent 
charge in free space has been determined. In the sisn,tlation the ratio E]^/E 2 
from the previous iteration Is used to calculate the effective charge density. 
The procedure is initiated by assuming |Ej^/E 2 | « 1 which yields the con- 

ventional result; the effective free space charge is 1/e^ the actual charge 
in the dielectric. 
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Fig. 4.4 Development of the equivalent charge in free spac 
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The charge distribution in Fig. 4.4(d) suggests that eq. 4,4 can be 
simplified to two surface integrals on planes located at y ■ ±d. 


-L 


W 


-W 


4itco 


/(x - x’)2 + (y - d)2 + (z - z')2 


/(x - x*)^ + (y + d)2 + (z - z')^ 


)dx’dz' 


I 

4irco 


W 


a(z') £n ~ + ~ (y d)2 + (Z - 2^)2 ^ 

-W (x - x') + *^(x - x*)2 + (y - d)2 -f (z - z*)i 


where W,L are the half width and half length of the material 


if (x ± L)^ » (y ± d)^ + (z - z’)^ we find 


4re( 


rw 

-w 


/ 0 / (y + d)^ 4 (z - 2*)^ N , 

) S,n _ d)2 + (z - z*)2 ^ 


assuming p(z’) ■ A + Bz* where A,B are constants 




ix'“l, 

dz' 

I X ' *-*L 


2(y - d) atan(~-j) - 2(y + d) atan(- -|- -^) ) - 


^ (((y - d)2 + s2)(in (s2 + (y - d)2) - 1 ) - 


((y + d)2 + s2)(in (s2 + (y + d)2) - 1))) 


s«2-w 


S"Z+W 


(4.7) 


The electric fields are related to the scalar potential by 


-V $ 
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Vfhich yields 


E„ - 


^ (2(A + Bz)(.t.n(-^) - .t.n(— 


- B((y + d) )ln((y + d)^ + s^) - (y - d) S.n((y - d)^- + s^))) 


(4.8) 

z-w 

g»z4v 


^ ((A + Bz) t jjj + ih 


+ 2B ((y + d) atan( y -" ~ ) - (y - d) atan(--^))) 


s*z-w 


' S“2+W 


(4.9) 


The actual charge distribution p cannot in general be represented in 
the form A + Bz. The fields can, however, be found by the superposition of 
solutions for each segment of a piecewise linear approximation to the actual 
charge distribution. The piecewise linear model Is chosen to be continuous 
and to have a mean square error less than a set fraction of the average 
charge density. Figure 4.5 shows an example of the segmented approximation 
chosen to satisfy a mean square error constraint and the original charge 
distribution it represents. The position of the segments are chosen in such 
a fashion as to shift the end points between successive Iterations thus 
minimizing the accumulative effect of discontinuities in the slope on the 
charging model. 

To summarize, the actual charge distribution is transformed using eq. 
4.6 to give an equivalent charge In free space which is then approximated by 
a continuous piecewise linear function. The total electric field is then 
found by summing the contributions of each segment calculated using eqs. 4.8 


and 4.9 
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4.1.2 Incident Current Profile 

As the surface potential reaches a level comparable to the beam 
accelerating voltage the electric field above the material strongly in- 
fluences the incident beam profile. An understanding of the beam dynamics 
can be obtained by tracing individual electron trajectories. 

The necessary equations of motion are; 


d^z/dt^ - -|e/m| E^Cz.y) 
d2y/dt2 - -je/roj Ey(z,y) 


(4.10) 


It is assumed that E is constant in time for the duration of the trajectory. 
This coupled system of equations was numerically solved using a third order 
Runge-Kutta method with adaptive time increment. The formulation is suited 
for second order differential equations for which a first derivative is not 
involved. The related difference equations for the z component are f Weher 
1967]: 


z(T + t) - z(T) + t(z(T) + 1/4 (kg® + k^h) 

z(T + t) - z(T) + 1/4 kg® + 3/4 kjj^ 

where z « dz/dt 

t » a time increment 
hg‘^*» -fe/m|t E (z(T), y(T)) 
kjj^- -|e/mlt Eg(z*, y*) 
z*- z(T) + 2/3 z(T) + 2/9 kg® t 

(a similar set of equations is used for the y-component) . 

The neglected terms are of the order t^. A measure of the local error 
can be obtained by introducing kg^ . 


k^^ - -|e/m[t E^(z(T + t), y(T + t)) 








58 


The error In z and z are given by: 

Error(z) - t(-l/4 + 1/4 

Error(z) - 1/2 k ° - 3/2 k ^ + k ^ 

e z z 

The method, incorporating an adaptive time increment procedure to maximize 
t eubject to local error estimate constraints on z,z,y,y, requires only 
two function evaluations per coordinate per iteration. A further check on 
the global accuracy of the trajectory subroutines was performed at charging 
equilibrium wherein the conservation of energy condition (eq. 3.9) was veri- 
fied across the surface. 

The trajectories are terminated on a plane defined by the upper surface 
of the material using an iterative selection of the final time Increment. 

By sampling the material's surface with test electrons we are able to 
deduce the incident beam profile. The Impact parameters, namely the angle 
of incidence and the kinetic energy of the Incident electrons as functions 
of position are obtained directly from the trajectory calculations. The 
incident current density Is obtained by comparing the relative deflection 
of neighbouring trajectories as indicated in Fig. 4.6. If we assume neigh- 
bouring electrons with the same initial velocity and a slight lateral dis- 
placement do not cross, then the net current through the upper and lower 
surface will be equal. Because there is no charge accumulation in the volume 
V, it follows that 

>Jyj^/jy2 " ^2^^1 (^*1^) 

Equation 4.11 is used to relate the incident current density to be used 
in eq. 4.2, to the beam current density via the calculated trajectories. 

4.1.3 Net Charge Accumulation 

The net current density, as indicated in eq. 4.2, represents the 



differtncc bctwtcn incoming and outgoing currtnta. Tha outgoing currants 
Includa tha backacatter currant, the secondary emission current and the bulk 

a 

conduction current. Although each has bean described in Chapter 2 the func- 
tional forms used in the simulation will be suimnarized hara. 


Jbi ■ 

O.i 

for Teflon and Kapton 

'^se 

1.55 exp (2(1 - cos0)) 

for Teflon 

a 

0.68 exp(2(l - cos6)) 

for Kapton 

Jc " 

3 

for Kapton 


■ 0 for Teflon 

where E - y component of electric field (V/cm) 

KE- incident kinetic energy (keV) 


In keeping with the discussion in Chapter 3, all of the results for Teflon 
are assumed to be secondary emission limited. 


4,2 Internal Charging Model 

The results presented in Chapter 3 indicate that the maximum internal 
fields at equilibrium occur at the front, face and in the non-lrradiated bulk 
of the material, The bulk fields can be obtained directly from the external 
charging model as they are outside the free charge distribution. The front 
face field can be estimated assuming a uniform dose rate throughout a volume 

defined by the mean electron range. The proposed charge distribution is 

* ■- . * ■ 

illustrated in Fig. 4.7. The net charge as used by the external model is 
~(pl P2^ where is a surface-depletion charge due to secondary emission 
, and p 2 Is a charge layer due to the primary electrons. The external surface 
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Fig. 4.7 Internal charging model for Teflon. 


A bulk conduction current to ground Is Included 
for Kapton. 
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e^E2 - Ej^ - pj/ C q 
E3 - E2 ■ 


(A. 12) 


Balancing the currents out of the first charge layer Pj^, we find 


dpi/dt - " ‘^se ■ 8^2 

dpi/dt - - gE^/cj) - (g/e^e^)Pj^ (4.13) 

Over the period of a time increment in the external charging model E^ 
will be approximately constant and eq. 4.13 reduces to an ordinary differen- 
tial equation with constant coefficients which can be easily solved. 


Pl<t) - (Jge " 

(4.14) 

(p^(0) - e^s^/g(Jgg - gEj/£^))exp(-g s/e^ep 

where T <, s <. T + t 

Pl(0) ■ charge at the beginning of the iteration 

The steady state solution is 

Pi - (Jgg - gEi/er)eoEr/g (4.15) 

^2 = Jgg/g - IEj^I (4.16) 

Note the field between the two charge layers E2 is independent of the external 
field and is equal to the expression obtained for the front face field in 
the one-dimensional analysis (eq. 3.27). 

4.3 Transverse Conduction Current 

The possibility of strong tangential fields near the sample edge suggests 
the inclusion of transverse currents and subsequent free charge transfer in 
the model. A three-dimensional analysis can be developed based on T^ig. 4.8 




D<t 


I 

h 

k 

i 

f. 

K 


I 


and eq. 4.1. If q is the charge contained in the volume In Fig. 4.8 

then 

-d,/dt - (J, - + (Jy - Jy’>Vo 

+ - J.’>Vo 

For an infinite strip J - J ' ■ 0. The change in the free charge per unit 

X X 

surface area is 

-dp/dt ■ (J - J ') + (j - J ’)y_/z (4.17) 

^ y y z zoo 

where p ■ q/x z 
0 0 

Assuming that the non -irradiated volume is blocking the transverse currents 
must be confined to a shallow surface layer with a thickness approximately 
equal to the electron range. Equation 4,17 can therefore be written 

8p/3t - Jy' - Jy + g7 (E^ - E^»)/z^ (4.18) 

where (E - E ’)/z • 9E /9z 

z z o z 

r « average electron range 

assuming 9E /9z > E /g • 9g/9z 

Z Z 

Tangential currents were not included in the general model as consider- 
able uncertainty persists concerning the transfer of charge from the material 
edge. It was felt that the edge could be characterized by two extremes: in- 

finite edge resistance and zero edge resistance. Infinite edge resistance 
is used to describe the boundary condition that *0 at the edge. Zero edge 
resistance is used to describe the boundary condition that 3J^/32 5 0 (i.e. 
no charge accumulation occurs at the edge due to surface tangential currents). 
Results have been obtained for both extremes. No attempt was made to model 
an inherent surface conductivity. 


i 
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5. SIMPLIFIED ANALYSIS 


The charging behaviour of the central portion of the dielectric strip 
when exposed to a normally incident uniform electron beam can be predicted 
using approximate field expressions. The approximations are valid when 
trajectory displacements from the centre are small compared to the half 
width of the material and are based on an assumed uniform charge distribution. 
The results are sufficiently accurate to be useful in practice and also 
provide data for later comparison with the fully numerical approach in 
Chapter 6. 

5.1 Central Fields of a Uniform Charge Distribution 

The following expressions are obtained from eqs. 4.7, 4.8 and 4.9 
assuming y » d and a uniform charge distribution (i.e. A ■ Pq, B • 0) 

a -— 2 -. (atan (^ ‘ ■" - atan( - - ) (5.1) 

V y 


.-isiav 

TTC« 


(g2 _ y2 » 


(y2 + (z - w)2)(y2 + (z + w)2) 


(5.2) 


Pod 


4w 


yz 


TTEo Cy2 + (z - w)2)(y2 + (z + w)2) 


(5.3) 


It is apparent from these equations that if the equilibrium condition is 
governed by the external fields, as is the case for emission-limited charging, 
an equilibrium state can be defined for an arbitrary material thickness by 
the steady state dipole moment per unit area p(z)d which produces a unique 
set of eixternal equilibrium fields. The steady state charge density is thus 
found to be Inversely proportional to the material thickness as implied in 
eq. 3.5 for a given equilibrium surface potential. 


Equations 5. I, 5.2 and 5.3 can be further simplified if we restrict our 

2 

attention to the region where z/w « 1. Ignoring terms of the order (z/w) 
we find 

^ • 2 p^d/ircQ atan w/y (5.4) 

Ey ■ pod2w/TTeo l/(y^ + w^) (5.5) 

(5.6) 

The equation of motion in the y direction, using eq. 5.5, is found to 

• 2 *2 

be Independent of the z coordinate. If we assume that (y) » (z) (i.e. a 

paraxial approximation) we find 

y “ ►'^ZBE/m - 2l$(y)e/m| (5,7) 


Equation 5.7 can be integrated numerically to find the vertical component of 
the trajectory as a function of time. The equation of motion in the z 
direction is more complicated but can again be treated numerically once y(t) 
is known. In general it takes the form 

z ■ f(t) z (5,8) 

where f(t) is a function of time defined 
by eqs. 5.6 and 5.7 

z( 0 ) » Zg , z( 0 ) - 0 

By simple substitution it is evident that if a solution to eq. 5.8 exists 
g]^(t) then a scalar multiple g 2 ■ <xgi(t) will also satisfy it with the new 
initial conditions: g 2 ( 0 ) “ agj(O), 82 ( 0 ) “ 0* Because the y component is 

Independent of z, two neighbouring trajectories represented by gj and g 2 are 
simply related by scaling the z coordinate of the path as indicated in Fig. 
5.1. The incident current density is therefore constant in the central 
region and can be determined by a single trajectory calculation as indicated 






in eq. 5.9 


, giCOH. - I) ii(0) 

•'A ■ ' tifo;. i ' > ■ <5-’> 

for any scalar a 

In addition since y is fixed at the surface ai^d z is found to be proportional 

I i 

to Zf the angle of intpact for small angles will also be proportional to z. 

Under symmetrical charging conditions the electrons incident at the 

centre of the sample will be directed normal to the surface. If the charging 

is emission-limited the central surface charge density can be calculated using 
» 

eqs. 5.3 and 3.13. 

Po - irco(BE - KE2)/(2dlel atan (w/d)) (5.10) 

where KE 2 " unity emission energy for normal Incidence 

For w » d eq. 5.10 reduces to the expressions obtained In Chapter 3 (eqs. 

3.5 and 3.13) describing a one-dimensional charging model. 

The trajectories for different beam energies are closely related. For 
a surface potential (and hence a charge density) equal to specific fraction 
of the beam accelerating voltage the resulting trajectories will be invariant. 
If y^ Zj describe a trajectory from a beam of energy BE^ onto a material with 
a surface potential equal to (k^BE^/e) then a corresponding trajectory y 2 » Z 2 
exists for a beam energy BE 2 and a surface potential (k^BE 2 /e) such that 
(from eq. 5.5): 

given (yj^ + w^)yj ^ 

( F 2 ^ + w ^)^2 - 

where kj * a constant 
y2(0) - yj(0) 
y2(0) - ^E2 /BEi ^^(0) 
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If we let 72 ■ yiCi/^E^TBE^ t) eq. 5.11 is satisfied along with the appro- 
priate initial conditions. From eq. 5.6 and using 72 * ** above, ^2 
found such that: 


given Zj/Zj - fj(t) 


Z2/Z2 - BE2 /BEj fj(/BE2/BEj t) (5.12) 

where Z2(0) - Zj^(O), ^(0) - Zj^(O) “ 0 

Substituting Z2 ■ z^(v'iE^7BE^ t) satisfies eq. 5.12. The resulting path 72, 

Z2 is unaltered as both coordinates are affected by an Identical time scaling. 
Using this result and calculating the incident current density as Indicated 
in eq. 5.9 we are able to generate a universal set of curves describing the 
incident current density as a function of the surface potential, as can be 
seen in Fig. 5.2. The accuracy of the analysis becomes uncertain when the 
surface potential approaches the accelerating beam voltage as the paraxial 
approximation used is no longer applicable. 

The angle of Incidence at equilibrium as a function of beam energy is 
given in Table 5.1 for Teflon. The equilibrium trajectories are more strongly 
deflected for higher beam energies as the equilibrium surface potential 
represents a larger fraction of the accelerating beam voltage. 

Table 5.2 presents the angle of incidence as a function of normalized 
coordinate z/w. The consistency of this value suggests a degree of width 
scaling is inherent in the equations of motion governed by eqs. 5.4, 5.5 and 
5.6. If the trajectory y^^, z^^ onto a sample of half width equal to unity 
is known, then it can be shown by direct substitution that 72, Z2 satisfy eqs. 
5.5 and 5.6 for an arbitrary halfwidth w when 


*2 " 

72 “ wyi(t/w) 


(5.13) 



Fig. 5.2 Reduction of Incident current density due to 


bd«m spreading 








with tht following initial conditions: 

22(0) ^ w 2 i(o) * " 

* 2 ( 0 ) - zj;(0) y 2 ^ 0 ) - yj(0) 

If the Dst trial is vary thin, we approximate the material's surface by the 
plane y ■ 0 which implies 


lyg-O 


(5.14) 


ly2-0 


At the material surface the energy of impact and hence the surface potential 
as well as the angle of Impact can therefore be represented as universal 
functions of z/w. 

The results in Fig. 5.2 were used to modify eq. 3.16 to account for the 
effect of beam spreading on the charging rates in Teflon. 


J^/J^^d - BS - SE) 


(5.15) 


where is given in Fig. 5.2 


The charging dynamics for beam energies from 5-25 keV based on eq. (5.15) 
and experimental data produced at Lewis Research Centre [ Purvis et al. 1979] 
for Teflon are compared in Fig. 5.3. NASCAP predictions and Lewis Research 
Centre data a« compared in Fig. 5.5 [ Roche and Purvis. 1979]. Both numerical 
simulations, as well as those reported by Beers . [1980] tend to be slower 
than the Lewis Research Centre data. Purvis et al [1977] have also reported 
excessive leakage currents in Teflon implying an effective conductivity of 
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Fig. 5,3 Charging dynamics for Teflon with beam spreading. 
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Fig. 5.4 Charging dynamics for Kapton with beam spreading. 
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•16 

10 (Hem) . The tuthore lift surface or edge leatuige as possible parallel 
conducCing paths to account for this observation. 

Charging results incorporating beam spreading are presented in Fig. 5.4 
for Kapton under conditions that support either eaission or conduction-limited 
equilibrium. No dramatic departure is observed in the general qualitative 
behaviour of either curve from those already discussed. 

Discharge studies [Balmain and Hirt . I960] on 50 pm thick Teflon irradiated 

2 

by 20 keV electrons over a range of current densities 0.5 - 50.0 nA/cm have 
indicated typical waiting times between discharges range from 60,0 - 2.5 
min. respectively. The corresponding range in the current denslty-tlme- 
thickness product is 9.0 - 37,5 nC/cm. Comparison with Fig. 5,3 suggests 
chat discharges occur after the net charge accumulation process has reached 
a steady state. 

5.2 Surface Fields 

The electric fields above the surface of a charged dielectric are cal- 
culated to determine whether such fields could initiate a breakdown, par- 
ticularly near the material edge, and to characterize the immediate environ- 
ment in order to determine its effect on the emission processes. The 
tangential field is continuous in the vertical direction whereas the normal 
electric field is discontinuous at the surface as noted In eq, 4.12. Both 
fields reach a maximum in the neighbourhood of the edge. 

For a uniform charge distribution the fields near the positive edge 
(z " w) and near the surface (y ■ d) can be calculated with the following 
definitions: 


4z 

m 

z 

- 

w 

2w 

a 

z 

+ 

w 

Ay 

m 

y 

- 

d 

2d 

m. 

y 

+ 

d 
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where ^2 ^ 

Ay « d , Ay i 0 
d « w 

The normal aurface field outside the charge distribution can be calculated 
from eq* 4.8 and is plotted in Fig. 5.6 versus the horizontal displacement 
from. the edge Az with the height above the surface Ay as a parameter. The 
position and magnitude of the maxima can be determined to be 

Az ■ f *^2dAy (5.17) 

Ey ■ pQ/2ireo(7r/2 ± d/w - »/2Ay/d ) (5.18) 

The field strength is well behaved in that no singularity Is encounter- 
ed as the surface is approached. The actual charge distribution is not pla- 
nar having a finite depth and an appropriate choice of Ay must be made to 
account for this. Fortunately, for sufficiently thick samples the maximum 
field is insensitive to Ay. 

If we consider the charge distribution in Fig. 5.7 the edge normal field 
using eq. 4.8 is 

E« • ((A + Bw)(d/w - atan Lz / atan Az/d) 

^ ^ (5.19) 

- Bd £n((d- + Az^)/w^)) l/2ireo 

where A - (p^ax + Pnin)/2 " Pavg 
B " <Pmin - Pmax^^w 

The normal electric field at the edge is reduced to zero when 

Pmin * Pavg^l ^1/(1 + 2£n(w/d))) (5.2Q) 

For pjjj^ less than this value the zero crossing moves in from the edge towards 
the centre of the strip and the normal fields at the edge are reversed tending 
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to attract ealsslon electrons back to the surface. For typical dimensions 

(w •• 0.375 cm, d ■ 0.01 cm) the zero crossing coincides with the material 

edge when p /p . - 1,28. If we Increase the ratio to p_ /p , - 3.0 
DAx nin uAx nin 

the surface normal field will be zero at ^ ■ -0.165 L (i.e. inside the 

z y 

edge) . Since represents the actual depth of the volume charge distri- 
bution (of the order of microns) the section of material under field re- 
versal is found to be very narrow for even large charge differentials. 

This coupled with strong edge tangential fields* to be described later in 
this section, suggest that the return of emission electrons to an Isolated 
sample should be an insignificant effect. 

For non-uniform specimens (where severe differential charging due to 
geometry or material composition could occur) the external fields can 
effectively ellmir.Ate the secondary emission and/or photo-emission process 
f Hand el et al . 1978] . This is particularly evident in the irradiation of 
strips of material possessing opposite charging characteristics (charging or 
non-charging) when placed alternately side by side. Figure 5.7 illustrates 
a subection of such an array where the non-charging strips form gaps that 
are assumed to be at zero potential. The equlpotential lines sketched above 
the surface suggest the development of a saddle point potential barrier 
representing a minimum energy requirement for an electron emitted from the 
gap to leave the surface. Since secondary emission and photo emission 
electrons tend to be low energy (” 10 ev) these processes can be surpressed 
by a relatively low saddle point potential. 

The potential above the centre of a grounded gap of width 2w between two 
uniformly charged samples each of width 2V? can be calculated using eq. 5.1: 

$(y) ■ -ppd/e^(atan(w/y) - atan(w + 2W) /y)/(tT/2) (5.21) 
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where PqA/Cq ■ surface pocenclal at the centre of the charged 
elements 

The saddle point is located by finding the root of the derivative of 

eq. 5*21 

yt - W/2w/W + w2/W2 (5.22) 

where y^, ■ the ' ^^ht aboyc the ground plane 

at which the saddle point is located 

Substituting eq. 5.22 into eq. 5.21 the saddle point potential barrier can 
be calculated. 

♦ ■ V(atan s - atan l/s)/(Tr/2) (5.23) 

where " potential barrier 
s ■ /I + 2W/w 

V ■ maximum surface potential 

Equations 5.22 and 5.23 are plotted in Fig. 5.8. As the fraction of the to- 
tal surface covered by the charging material is decreased the potential 
barrier is reduced in magnitude and moves away from the surface. Excellent 
agreement has been obtained between the above analysis and the NASCAF program 
results [Stannard et al . 1981]. 

The tangential field near the edge of a uniformly charged strip can be 
calculated using eq. 4.9 and has the following form: 

<a/w)2)) (5.24) 

I^ote as both A2 and Ay approach zero the tangential surface field has a 
logarithmic singularity. The maximum field is located directly above the edge 
(Az " 0) and is given by 

Eg - Po/27reo ()in(?;d/Ay) - 1/2 inXl + (d/w)S) 


(5.25) 
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If we assume that Instead of a planar charge layer there exists a uniform 
volume charge distribution of depth r, then the maximum tangential field can 
be obtained by Integrating eq. 5.25 with respect to Ay. The resultant 
edge field at the centre of the charge distribution Is 

Eg - Po/2tt6o (An Ad/r + 1 - 1/2 An(l + (d/W)^)) (5.26) 

For sufficiently thick saiftples eqs. 5.26 and 5.25 are the same If we choose 
Ay • r/5,5, Since the details of the Internal charge distribution are not 
known we have adopted the convention that all "surface" fields will be evalu- 
ated at Ay « 0.01 d. Rescaling of our results for an alternative charge 
distribution can be accomplished uslnf eq. 5^25. 
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6. TWO-DIMENSIONAL NUMERICAL RESULTS AND DISCUSSION 


The numerical results have been chosen to lllusi:rate the charging 
response to a variety of geometries and beam characteristics. The Teflon 
simulations were assumed to oe emission-limited and therefore Independent 
of beam current density. A standard current density of SO nA/cm was 
selected for these runs. Data was accumulated at 160 points across the 
surface. Trajectory calculatloi^s were made with a nearest-neighbour 
spacing ranging between 0.25% and 5.0% of the total sample vrldth, typically 
requiring 50 trajectories per iteration, vd.th a maximum resolution being 
obtained in the neighbourhood of the material edges. Standard runs 
required 8 min. of CPU time on the University of Toronto Computing Services' 
IBM 3033. The piecewise linear charge approximation was chosen to have a 
root mean square error less than 1.3% of the average charge density. 

Under symmetrical charging conditions the expected solution symmetry 
was not Incorporated into the charging model. Instead the full surface 
was represented so that asymmetric discrepancies can be used as a measure 
of the model’s reliability. 

The net charge density refers to the net negative charge in the sample. 
The surface charge density refers to the accumulation of an assumed 
positive charge layer on the material surface due to secondary emission. 

The surface fields refer to fields calculated a distance 0.01 d above an 
assumed planar (zero-thickness) distribution of net charge. 

6-1 Beam Energy Dependence for Teflon 

The equilibrium net charge profiles for Teflon irradiated with beam 
energies ranging from 5 to 25 keV in 5 keV steps are given in Fig. 6.1 
(material half width 0^75 cm, material thickness 100 imi). In Fig, 6,2 
the equivalent charge in free space as described by eq. 4.6 is presented. 

In general, the internal normal bulk field is much larger than the 
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Fig. 6.1 Equilibrium surface charge density distributions for Teflon. 
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Fig. 6.2 Equilibrium equivalent surface charge distributions 
in free space and piecewise linear approximations 
used in field calculations. 
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external normal surface field and the equivalent charge is therefore 
approximately equal to 1/e^ of the actual charge. The scales in Figs. 6.1 
and 6.2 differ by a factor of 1/e^ in order to facllatate this comparison. 
The crosses in Fig. 6.2 represent the end points of the linear segments 
used to approximate the effective charge density for the purpose of 
making field calculations. 

The charging dynamics are illustrated by comparing intermediate stages 
in the evolution of the net charge in Figs. 6.3, 6.4 and 6.5 for beam 
energies of 20, 10 and 5 keV respectively. The vertical scales in these 
figures were chosen to be nearly inversely proportional to the beam energy. 
To a first approximation the curves for beam energies greater than 10 keV 
have the same shape indicating that a rather simple beam energy scaling of 
the equilibrium charge distribution is applicable to the entire surface. 

A more quantitative analysis in Fig. 6.6 suggests that the equilibrium 
charge density across the sample is proportional to (BE-k) where k is a 
function of position on the surface. At the centre k is equal to the 
second total emission crossover energy at normal incidence (2.1 keV for 
Teflon) as suggested by eq. 3.13. The value of k is approximately 
constant over central portion of the surface and reaches a minimum of 
1.05 keV near the edges. 

The charging dynamics at the edges differ considerably from those at 
the centre where the net charge density monotonlcally approaches a steady 
state in a manner consistent with eq. 5.15. At the edges the charge 
density overshoots its equilibrium value. This is illustrated in Fig. 6.7 
where the temporal behaviour of the charge density at the centre and edges 
is explicitly presented for a 20 keV charging beam. Comparison with the 
predicted centre response from Fig. 5.3 indicates excellent agreement. 

The overshoot at the edges tends to increase with increasing beam energy. 
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Fig. 6.5 Charging dynamics £or,Teflon. Beam energy 5 keV. 
Beam density 50 nA/cm“. 
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Fig. 6.7 Explicit time response at centre and edge of sample. 
Beam energy 20 keV. 
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A possible sxplanscion for this bshavlour is prasantad in a discussion 
of the nactssary squilibrium conditions later in this section. 

The steady state surface potential distribution for each bean energy 
is presented in Fig. 6.8. At higher bean energies the equilibrium surface 
potential becones less uniform as the surface charges to a larger fraction 
of the bean accelerating voltage and the beam is more strongly deflected. 
Surface potentials measured by Robinson [1077] and Stevens et al [1978] 
are also presented in Fig. 6.8. Considering the uncertainty in the experi- 
mental points as well as the differences between the physical and model 
geometry the agreement is quite good. 

The angle of impact in Fig. 6.9| as suggested in Chapter 5p displays 
a near linear central region, particularly for the lower beam energies, 
which increases in slope for increasing beam energies. In general the 
slope is larger than that found in Table 5.1; the reason may be that the 
non-uniform charge distribution tends to increase the tangential fields 
and enhance the degree of beam spreading. 

The angle of impact is closely associated with the surface potential 
through the current balance in eq. 3.10 and the conservation of energy 
expressed in eq, 3.9 

KE + lev] - BE 

for emission-limited equilibrium 

J^(l - BS - SE) - Jc * 0 

therefore 

BS(KE,0) + SEOCE.0) - 1 (6.1) 

where“KE,0 atvl V are unknown 

An explicit relationship between the equilibrium kinetic energy 
and the angle Of Impact is given in Fig. 6.10 for a 
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Fig, 6.8 Equilibrium surface potential distributions for Teflon, 

>f experimental data points from Stevens et al , [1978] 
(experiment: Teflon strips 0,013 cm thick>halfwidth 
S . 0 cm) 

0 experimental data points from Robinson , [1977] 

(experiment: circular Teflon samples 6.013 cm thick 
defined by an aperture 2.5 cm in diameter cut in a 
1.3 mm aluminum plate placed in contact above the 
material) 







Fig. 6.9 Angle of incidence at equilibrium. 



Fig. 6.10 Unity emission energy vs. angle of incidence at 
equilibrium for Teflon. 
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total emission coefficient as described by Wall et al [1977] (eqs. 2.1, 

2.2, 2.7 and 2.9) and as described by Robinson and Budd [1980] (eqs. 2.8 
and 2.10). Had the angular dependence suggested by Robinson (eq. 2.10) 
been used in the computer simulation the central region of the sample would 
have been more uniform. Figs. 6.8 and 6.9 compare favourably with the 
predicted relationship in Fig. 6.10. 

The charging dynamics can be explained in «*:erms of the relationship 
between the surface potential and the angle of Impact. As the sample begins 
to charge the beam is largely undeflected and a uniform charge distribution 
is developed. However, as the surface potential approaches the accelerating 
beam voltage the angle of impact at the edges reaches a level (approximately 
50°) where the associated local equilibrium potential begins to drop 
signifigantly, forcing a loss of charge primarily through an increase in the 
secondary emission coefficient from the edge. This loss of charge enhances 
the tangential fields thereby maintaining a large angle of Incidence. As 
a result the charge density at the edge collapses as the potential at the 
centre nears the beam Voltage. Since a 5 keV beam charges Teflon to a 
maximum surface potential of approximately 3 kV (representing only 60% of 
the accelerating voltage) an overshoot does not develop and the resulting 
charge and surface potential distributions are more uniform. 

The external normal surface field, external tangential surface field 

and the internal normal bulk field distributions are given in Figs. 6.11, 

6.12 and 6,13 respectively. The maxima obtained by the external normal 

surface fields are in general less than that required for dielectric 

breakdown. The edge surface tangential field and the central internal 

6 

normal field are comparable both exceeding 10 V/cm for the 20 keV and 
25 keV beams. For a strip width of 1.5 cm the tangential field maxima 
are 60% smaller than the peak internal bulk field. 
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Fig. 6.11 External normal surface field distributions at equilibrium 
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Fig. 6.12 External tangential surface field distributions at equilibrium 
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The strong peak In the tangential field has a half-maximua width of 
the order of the sample's thickness, As indicated by eq. 5,25 the peak 
tangential field is influenced by the assumed internal volume charge 
distribution and is therefore subject to some uncertainty. In many 
experimental setups the sample is also covered by a grounded metal aperture 
plate which could enhance the peak surface tangential component at the edge. 

The tangential field calculated at y - d/2 (halfway through the material) 
is approximately 1/5 the magnitude of the surface tangential component. The 
undulations in the surface fields are due to the discontinuities In the 
slope of the piecewise linear charge representation. The actual quantities 
plotted in Figs. 6.11, 6.12 and 6,13 are the magnitude of the various field 
components. The direction of the tangential field should be understood to 
be towards the centre of the strip and passes through a zero at that point. 

The paths traced by low energy secondaries emitted from the surface 
charged to equilibrium by a 20 keV beam are approximated in Fig. 6.14. 

The electrons are assumed to be emitted normal to the surface with an 
energy of 1.0 eV. A small fraction (approximately 7%) are found to be 
turned back to the ground plane. Experiments performed by Cuchanski [1978] , 
in which an electron microscope was used to image the secondary emission 
occuring from a charged dielectric, indicated a band near the grounded 
metal edge from which no secondaries reached the detector. The band 
measured approximately 9% of the sample radius in good agreement with the 
predicted behaviour. The blowoff electrons during the discharge event 
cannot be inferred from Fig. 6.14 due to severe space charge limiting that 
has been ignored in this treatment f Stettner, 1980; Katz, 1980] . Measure- 
ments made by Balmain and Hirt [1980] indicate that the fraction of the 
released charge returning to the ground plane is increased to approximately 
50% during the discharge event. 
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Figures 6,15, 6.16 ind 6.17 Indicate the incident current denaity 
profile, the dose rate profile and incident trajectories taken at steady 
state for the sanple charged by a 20 keV bean. The nornalized current 
density has a broad minimum of 0,53 at the centre of the strip and maxima 
at the edges of 0.7. The central current density is in good agreement 
with the simplified analysis presented in Fig. 5.2. The dose race is 
remarkably flat considering Che variation in its constitutive parameters: 
angle of incidence, surface potential and current density^ as given by the 
formula 

^ » / r, .>-.-0,358 

D a Jj^/cosQ KE 

which is derived from eqs. 2.12, 2,13 and 2,14. 

The internal front face field is presented in Fig, 6.18. The equili- 
brium field is governed by the relationship (eq, 3.26) 


'ff 


-3 /g 
se ® 


It was demonstrated (eq. 3.27) that this quantity is a rather weak function 
of incident current density and impact energy so the uniformity of the two- 
dimensional results should not be too surprising. The magnitude of the 
field agrees quite well with that predicted in Chapter 3 when the degree of 
beam spreading is taken into account. 

[Effl- (J^)°*^ 6.1 X 10^ V/cm 

■ 3^ • 0.53 where 3^ ■ 50 nA/cm“ 
therefore the one-dimensional prediction is 

- 1.64 X 10^ V/cm 


The surface charge layer required by Gauss’ Law to match the internal front 
face field and the external normal surface field in Fig, 6,18 and Fig. 6.11 








• lUf unwmi&'i! 


NORMflLIZED POSITION (X/HI 

Fig, 6.15 Reduction in incident current density due to beaiti 
spreading at equilibrium, iBeam energy 20 KeV, 
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Fig, 6,16 (Dose rate/beam current density) distribution at 
equilibrium for Teflon charged by a 20 keV beam. 
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Fig. 6,18 Intenial front face field for Teflon at equilibrium, 
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Fig, 6.19 Front face charge density at equilibrium 
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respectively is presented in Fig. 6.19. In generel e snell positive 

surface charge is required becoming negative at the edge due to the obeerved 

aaxioa in the external normal surface field* 

The average charge density at steady state is plotted against the 

equilibrium surface potential at the centre of the strip in Fig. 6.20. 

The upper curve is obtained from the one-dimensional model which produced 

a capacitance per unit area of The lower curve is based on the 

two-dimensional charging model results. The temporal evolution of the 

average charge density is plotted against the central surface potential in 

Fig. 6.21a. The charga-voltage relationship is found to follow the one- 

dimensional curve until near steady state when charge is lost from the 

edges. Comparable experimental data [Stevens et al , 1977] is reproduced 

in Fig. 6.21b, 11 the low voltage portion of Stevens' data is extrapolated 

2 

to the steady state voltage we find (for an area of 300 cm , d ■ 127 urn and 
- 2.25) good agreement with the one-dimensional curve in Fig. 6.20. 

The experimental results hewever, indicate an equilibrium average charge 
density that is lower than the computer simulation results by a factor of 
1/1.3. Such a discrepancy could be attributed to a change in the dielectric 
permittivity due to the high internal fields, Interestingly if we take data 
for the charging dynamics of Teflon r Purvis et al , 1977] in Fig. 5. A and scale 
the time axis by a factor of 1.5 as would be appropriate for the proposed 
change in we again obtain a much improved agreement with the numerical 
results as indicated in Fig. 6.22. Note that at a lower beam energy 
(hence a lower internal field strength) reasonable agreement already exists in 
the charging dynamics (Fig. 5.4). 

6,2 Ilaterial Width Dependence for Teflon 

The equillbriiun net charge density profiles in Fig. 6.23 were obtained 
for Teflon with half widths W equal to 1,5, 1,5/2, 1.5/4 and 1.5/8 


ont-diaensional analysis 
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Fig. 6.20 Average charge density at equilibrium vs. central surface 
potential. Assuming ■ 2.25 
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Fig. 6.21 Charging dynamics for Teflon. 





103 


(■aterlal thickness 100 un, bean energy 20 keV), In the centre the charge 
distribution is nearly independent of the sample half width when plotted 
against .the normalized co-ordinate z/W, At the edges the narrower strips 
collect signiflgantly greater charge density. 

The surface potential distribution when normalized to the width of 
the sample is also found to be insensitive to the sample width as indi- 
cated in Fig. 6.24a . The close agreement between experimental results 
and the simulation results in Fig. 6.8 where the range in half width was 
0.75 to 5.0 cm tends to support this form of width scaling. Further 
experimental evidence presented in Fig. 6.24b confirms the 
scaling on strips of Teflon thermal control tape [Aron and Staskus . 19791 .The 
behaviour of the surface charge distribution is also consistent with 
measurements indicating that the total charge released during a discharge 
event is proportional to the exposed area [ Balmain and Dubois . 1979] . 

Equation 6.2 is an approximate empirical expression fitting the 
behaviour of the surface potential distribution for arbitrary beam enerf,y 
and arbitrary sample width. 


V = Vq (1 - exp(- '*/(4.6)‘+ (1.185 - Iz/wD** - k(W)» 


where 


( 6 . 2 ) 


Iz/Wj S 1.0 

k(W) ■ a correction factor to match the edge characteristic 


Vq “ equilibrium surface potential from eq. 3.13 
from Fig. 6.24 we find 

k(W) - 0.5120, 0.4945, 0,4464, 0.350 
W » 1.5, 1.5/2, 1.5/4, 1.5/8 cm respectively. 


for 


NORMflUZED POSITION (K/H) 


Fig. 6.23 Equilibrium surface charge density distributions fpr 
Teflon. 










For | 2 :/w| < 0,9 eq. 6,2 can be slopllfied to a function of a aingle 
nomallzed parameter z/W 

V a Vq (1 - exp (-4,6(1, 185 - |z/w|))) (6.3) 

where 

Iz/W| < 0,9 

Similar width scaling Is found in comparing the angle of Incidence as 
a function of the normalized coordinate for each sample width. This result 
could have been deduced directly from the equilibrium surface potential 
curves since the angle of Incidence at equilibrium Is fixed once the 
surface potential and beam energy are specified (using Fig, 6.10), The 
angle of Incidence for each scrip is therefore practically the same as that 
presented In Fig. 6.9 for a 20 keV beam. Slight shifts between the curves 
at Che edges were found in accordance with the noted divergence in the 
surface potentials in Fig. 6.24. 

The bulk field, being approximately equal to the surface potential 
divided by Che material thickness, is also found to obey the same width 
scaling. A representative curve can be found in Fig. 6.13 for the 
20 keV beam. 

The external surface fields are presented In Figs. 6.25 and 6.26. 

Near the centre of the samples they are found to be inversely proportional 
to Che strip width. The front face fields are Che same as those shown in 
Fig, 6.18. 

6.3 Beam Angle Dependence for Teflon 

The equilibrium net charge profiles in Fig, 6,27 were obtained for 
Teflon by charging with electron beamis incident at 0“, 15®, 30® or 45® 
relative to the surface normal (material half width 0.75 cm, material 
thickness 100 um, beam energy 20 keV) . The beam originates from Che left 
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of Che fpecinen tending to deposit • larger charge density at the edge 
nearest to its source. The central aaxiaum tends to shift towards the 
near edge and the maxinum charge density tends to increase with increasing 
beam angle. 

The charging dynamics for a !>eain angle of 45* are Illustrated in 

Fig. 6,28. The explicit time dependence for the centre and edges is given 

2 

in Fig. 6.29. Note the beam current density of SO nA/cm refers to the 
current component in the y direction (i>e. normal to the exposed surface 
not perpendicular to the beams direction). At the near edge the charge 
density had not reached an equilibrium value when the simulation was 
terminated as the rate of charge accumulation indicated by the final slope 
in Fig. 6,29 was still klL of the beam current. 

The surface potential distribution is presented in Fig. 6.30. The 
maximum is found to move to the near edge as the beam angle is Increased. 
The magnitude of the maximum remains constant and indicates the position 
on the sample where the incident electrons impact the surface at 0*. The 
angle of incidence is presented in Fig. 6,31 for each case. At the near 
edge 30^/3 z increases with increasing beam angle. For a beam angle of 45* 
30^/3 2 reaches 117“/nBn compared to an average of 8* /mm for a beam angle of 
0 *. 


The maximum internal bulk field as suggested by the surface potential 
does not change in magnitude though it moVes towards the near edge for 
increasing beam angles. Over a large fraction of the sample the bulk 
field is reduced. 

The external surface electric fields (Figs. 6.32 and 6.33) are found 
to be enhanced at the near edge and reduced at the far edge. The maximum 
tangential component was increased by a factor of 3 between the 0* and 
45* charging beams to 3.0 x 10^ V/cm. A similar Increase in the maximum 
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Fig. 6.32 External normal surface field distributions at equilibrium. 



Fig. 6.33 External tangential surface field distributions at equilibrium 
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Fig, 6.35 (Dose rate/beam current density) distribution at 
equilibrium. Beam angle 45°. 
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noraal cooponent resulted In a field strength of 6 x 10^ V/m for a 45'* 
charging beam. 

The front face field is found to be remarkably flat with no discenii- 
ble difference between the 45* and 0* beams, 

The incident current density and dose rate profiles at equilibrium 
for a beam angle of 45* are presented in Figs. 6,34 and 6,35 respectively. 

6.4 Material Thickness Dependence for Teflon 

The equilibrium net charge profiles in Fig, 6,36 were obtained for 
Teflon with material thicknesses of 100, 50 and 25 urn (material half width 
0.75 cm, beam energy 15 keV), The dipole moment per unit area p(z)d, as 
suggested in section 5.1, is independent of the thickness. The surface 
potential and angle of incidence are also Invariant. 

The tangential surface field (Fig. 6,39) is increased at the edge for 
a decrease in the material thickness. The external normal surface field 
(Fig, 6.38) shows no appreciable change. The internal bulk field (Fig. 

6,37) is Inversely proportional to the material thickness as would be 
expected for a constant surface potential. 

The front face field displayed no appreciable change from that presented 
in Fig, 6.18. 

6.5 N on-conventional Charging Geometries for Teflon 
6.5.1 Grounded Central Slot 

The effect of a grounded gap between two adjacent strips of 
Teflon has been investigated for gap widths of 0,283 mm and 0.660 mm. 

The geometry modelled is presented in Fig, 6.40. (material half width 
0.75 cm, material thickness 100cm, beam energy 20 keV) . Except for the 
peaks in the charge density at the inner edges the equilibrium charge 
distributions in Fig. 6.41 are similar to those found for a uniform 
material. 





Fig. 6.36 Equilibrium surface charge density distributions. 





Fig. 6.38 External normal surface field distributions at equilibrium. 
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Fig. 6.39 External tangential surface field distributions at'eqdiribftam 
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The surface potential distribution is presented in Fig. 6.42. 
The angle of impact (Fig. 6.43) suggests considerable beam steering in 
the neighbourhood of the slot. Figure 6.44 traces sample incident electron 
trajectories at steady state for the 0.66 mm slot width. The equili- 
brium current density is given in Fig. 6,45 indicating considerable 
beam focusing into the slot area, The dose rate profile at equilibrium 
is very similar to that obtained for a uniform strip. 

The external normal and tangential surface fields are presented 
in Figs. 6.46 and 6,47. Both display additional peaks at the inner edges. 
The narrower slot tends to produce a slightly larger maxima for both field 
components. 

The internal bulk field can be determined from the surface 
potential. Its maximum value Is little changed from that obtained for a 
uniform strip. The internal front face field is also similar to that 
found for a uniform strip. 

Calculation of the equipotential lines at steady state reveals 
a saddle point 0.9 mm above the centre of the 0.283 mm gap. The corres- 
ponding barrier potential was 13.2 kV. These results compare favourably 
with the value obtained from the approximate analysis of Chapter 5 
(eqs. 5.22 and 5.23 yield a saddle point potential 14.9 kV, 1,04 mm above 
the gap centre). 

Measurements made by Robinson [1977] suggest a large increase 

in the rate of breakdown occurs when a slot is cut in the material. A 

2 

300 fold Increase was observed for a 21 keV beam at 40 nA/cm when a slot 
was cut in the sample. The simulation results indicate that the Increased 
tendency to discharge is consistent with Increases in the surface fields 
and exposed edge length rather than the Internal fields. 









121 


6.5.2 Precharging Effects 

The effect of a non-zero Initial charge distribution on Teflon 

has been investigated for a uniformljr charged central sub-strip (material 

half width 0.75 cm, material thickness 100 urn, beam energy 20 keV) . A 

2 

negative charge region of 250 nC/cm was used at two different widths, 

0.94 mm and 0,235 mm. A positively charged region was investigated for 
the same charge density magnitude at a width of 0.47 mm. The charging 
dynamics are presented in Figs. 6.48, 6.49 and 6,50 for each case described 
above . 

The charge density distributions tend to converge to a universal 
steady state typified by those simulations with zero initial conditions. 
However, near the boundaries of the precharged region the charge distri- 
bution continues to reflect the discontinuous nature of the assumed 
initial state. 

The surface fields display the greatest sensitivity to the 
irregularities in the charge distribution at the end of the simulation 
(the term equilibrium is avoided as the net rate of charge transfer to 
the boundary region at the end Of the simulation was at least an order of 
magnitude greater than that found over the remainder of the surface) . 

The surface fields for the 0.94 mm negatively precharged stripe (Figs. 

6.51 and 6.52) are typical of the other cases described. The magnitude of 
the field pertubation is of the order of 10^ V/ct in both the normal and 
tangential component. 

A non-zero initial charge distribution can be physically 
attributed to a partial discharge of the material’s surface or static 
charging of the surface due to brushing or rubbing in the material's 
preparation. The charging dynamics in Figs. 6.48, 6.49 and 6.50 suggest 
that although the final charge distribution will approach that previously 




Fig. 6.50 Charging d^mamics. Strip 0.047 cm wide, precharged to 

-250 nC/cm^ (electron charge density). Beam density 50 nA/cm 
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Fig. 6.51 External normal surface field distribution above the 
final charge profile in Fig. 6.48. 
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Fig. 6.52 External tangential surface field distribution above the 
final charge profile in Fig. 6.49. 
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observed for zero initial conditions, the charging rates across the sample 
are characteristic of each region’s level of development. Near the 
precharge boundary the net charging rate approximates an average of the 
neighbouring rates causing highly localized perturbations that persist 
after the remainder of the surface has obtained equilibrium. 

Laboratory experiments in which dielectric films have been 
brushed before irradiation indicate that the direction of brushing can 
influence the behaviour of the visible discharge arcs [ K.G. Balmain , 

J.V. Staskus . personal communications] . Multiple arcs that parallel the 
direction of brushing were produced with either a camel hair brush or a 
paper tissue along either straight or circular paths. 

A possible mechanism for such behaviour can be developed 
based upon the characteristics described in this section. If the sample 
were charged along lines parallel to the direction of brushing then 
irregularities in the developing charge profile could enhance local 
external surface fields, thus Increasing the probability of discharge 
along these lines. If the discharge event is confined, to these high 
field regions subsequent charging would again develop such conditions 
along the edge of the discharge region thus perpetuating the area's 
discharge activity. 

6.5.3 Non-uniform Material Thickness 

The equilibrium net charge density profiles for a composite 
material consisting of two adjacent uniform sections i</ith sample thick- 
nesses of 50 urn wd 100 urn are indicated in Fig. 6.53 (sections full 
width 0.75 cm, beam energy 15 keV) , The general distribution is in 
agreement with those conclusions reached in determining the thickness 
dependence of the charging characteristics (section 6.3). 






The rctuXting surface potential (Fig, 6,54) however resembles 
that obtained when a gap was present In the material. The external 
surface fields (Figs, 6,55 and 6.56) also reflect the discontinuity of the 
material. Doth components have maxima ocurring at the centre and edges. 

The internal bulk field is presented in Fig. 6.57. 

6,6 Beam Current Dependence for Kapton 

The equilibrium surface potential for Kapton was found to be sensitive 

2 

to typical beam current densities (0.1 - 100.0 nA/cm ) for material 
thicknesses of the order of 50 urn and beam energies near 20 keV. As 
indicated in Fig. 3.2 there exist two distinct regions of equilibrium: 
emission-limited and conduction-limited. The emission-limited surface 
potential is independent of the beam current density and the material 
thickness as is the case for Teflon. Gonversely. the conduction-limited 
process is strongly influenced by both the current density and material 
thickness. These charging parameters have been chosen in this study to 
characterize the two-dimensional nature of each limiting mechanism. 

Simulations were conducted for thicknesses of 100 urn and 50 um. and 

2 2 

at current densities of 100 nA/cm and 1 nA/cm (material half width 0.75 cm 
beam energy 20 keV). 

The net charge density profiles at equilibrium are presented in 
Figs. 6.58 and 6.59. The high current thick sample simulation typifies 
the emission-limited state obtained with the Teflon model. The low current 
thin sample simulation is conduction-limited and displays a comparatively 
flat, uniform charge distribution. 

The equilibri'-an surface potential is presented in Figs. 6,60 and 6.61. 
The conduction-limited potential is constant while at the edges the 
emission-limited value decreases to less than 50 % of its peak central 
magnitude. As suggested in Chapter 3, Kapton if emission-limited tends to 
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Fig. 6,58 Equilibrium charge density distributions for Kapton. 
Material thickness 100 yro. 
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Fig. 6.59 Equilibrium charge density distributions for Kapton. 
Material thickness 50 pm. 
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charge to within 1 kV of the beam accelerating voltage compared to 2 kV for 
Teflon. The deflection of the beam at steady state is therefore greater for 
Kapton. The angle of incidence graphs (Figs. 6.62 and 6.63) indicate 
maxima at the edge of almost 80” compared to 65” for comparable Teflon 
simulations. 

The uniformity of the results for conduction-limited charging can be 
understood in terms of the simplified analysis presented Chapter 5. The 
incident current density onto a uniform charge distribution was found to be 
constant over the central portion of the strip as would be the bulk internal 
fields. Tt\e conduction current density, the main charge balancing mechanism, 
would therefore also be uniform and a consistent equilibrium could be reached 
under the assumed charge profile. 

The external surface normal fields (Figs. 0.64 and 6,65) are comparable 
to those obtained in the Teflon simulations. The tangential surface fields 
(Figs. 6,66 and 6,67) suggest that the flat conduction -limited charge 

9 

distribution produces slightly larger maxima at the edges than does the 

emission-limited charge distribution. 

The internal bulk fields (Figs. 6.68 and 6,69) have broad maxima at 

the centre of the strip. The internal front face fields (Figs, 6,70 and 

6,71) are generally constant across the surface. The incident current 

0 3 

density dependence suggested by eq, 3.27, ’ , is clearly evident. The 

required surface charge layers are presented in Figs, 6.72 and 6.73. 

6.7 Tangential Current Results for Kapton 

A tangential transfer of charge through an Irradiated front surface 
layer (as outlined in section 4,3) has been implemented for Kapton in the 
following charging simulations. Two separate edge boundary conditions 
were investigated; 3J /3x ■ 0 and J -0, The term zero edge resistance 

Z Z' 

is used to describe the first case where there can be no charging at the 



Fig, 6,62 Angle of incidence at equilibrium for Kapton. 
Material thickness 100 pm. 
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Fig. 6.63 Angle of incidence at equilibrium for Kapton 
Material thickness 50 ym. 
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Fig. 6.64 Normal surface field at equilibrium for Kapton 
Material thickness 100 pm. 
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Fig. 6.65 Normal surface field at equilibrium for Kapton. 
Material thickness 50 pm. 
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Fig. 6.68 Internal normal bulk field at equilibrium for Kapton. 
Material thickness 100 urn. 
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Fig. 6.73 Front face charge density distributions at equilibrium 
for Kapton. Material thickness 50 um. 
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edge due to transverse currents, The tern Infinite edge resistance is used 
to refer to the latter condition, The field gradient at the edge to be 
used in the calculation (cq, 4.18) is made between the edge and the 
adjacent material data point. 

The material half width and the beam energy were kept constant at 

0.75 and 20 keV respectively as was the case in section 6,5, Each boundary 

condition was tested on two models corresponding to emission-limited and 

conduction-limited conditions (material thicknesses 100 urn and 50 um, beam 

2 2 

current density 100 nA/cm and 1 nA/cm respectively) , 

The charging dynamics are presented at equal time intervals in Figs. 
6.74, 6.75, 6,76, 6,77, 6.78 and 6, 79. The results for an infinite edge 
resistance show no noticeable deviation from the curves generated for no 
transverse currents. With the exception of a slight depression at the 
edge the results for zero edge resistance also indicate no deviation. 
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Fig. 6.74 Charging dynamics for Kapton. Material thickness 100 ym, 
current density 100 nA/cm^, No tangential currents. 
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Fig. 6,7S Charging dynamics for Kapton. Material 
current density 1 nA/cm^. No tangential 


m. Material thickness 50 ym, 
lo tangential current. 
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Fig. 6.76 Charging dynamics for Kapton. Material thickness 100 pm, 
current density 100 nA/cm^. Tangential currents with 
boundary condition J^^O. 
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Fig. 6.77 Charging dynamics for Kapton. Material thickness 50 urn, 
current density 1 nA/cnr. Tangential currents with 


boundaiy condition J_*0. 
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Pig. 6.78 Charging dynainics for Kapton. Material thickness 100 um 
current density 100 nA/cm^. Tangential currents with 
boundary condition 3J2/3z*0. 
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Fig. 6.79 Charging dynamics for Kapton. Material thickness SO Mm, 
current density 1 nk/cvr. Tangential currents with 
boundary condition 572/32*0. 
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7. CONCLUSIONS 

A tvo-dimtnfion«l nodtl has baan prasantad to dascciba tha accuaula- 
tlon of charge on long, planar dlalactrlc strips by a unlfom nonoanergetic 
alactron baani used to simulata spacecraft charging conditions. Nunerlcal 
results detail the net surface charge density, surface potential and 
material field profiles for both Teflon and Kapton. Teflon has been 
treated as being completely emission-limited whereas Kapton is governed by 
both emission and conduction currents, The response of both materials under 
a variety of charging conditions has been investigated, 

The equilibrium state can be categorized as being either edge-like or 
centre-like. The edge, under normal charging conditions, is characterized 
by maxima in both the normal and tangential components of the external 
surface electric field and by minima in the surface charge density and 
surface potential. In general, the total external surface field at the 
edge is dominated by the tangential component. The centre is characterized 
by a maximum in both the internal bulk electric field and the net surface 
charge density, and by a minimum in both components of the external 
surface fields. 

The one-dimensional analysis of Chapter 3 adequately describes the 
behaviour of the centre region whereas a complete two dimensional model is 
necessary at the edge. In the course of this study a positive correlation 
has been established between the edge»llke * glons and the apparent sites 
of discharge initiation and blow off whereas punchthrough discharges can 
be associated with the high internal bulk fields at the centre. 

For emission limited charging and normal beam incidence it was found 
that the maximum surface and bulk fields were approximately proportional 
to the beam energy. The internal normal bulk field was found to be 
inversely proportional to the material thickness d and Independent of the 


macerlal halfvldeh U. Tht mcximum tdge tangtntial fltld wat found Co vacy 
—0 61 «0 33 

as (d) and (U) , Tht location of eha naximun fltld in cht naCtrial 

chertfoce dtptnda on Cht ratio W/d, If W/d 4 20 tht maxinun nattcial fitld 
it the surface tangential component at the edge; otherwise the maximum is 
Che internal normal bulk component located at the centre of the strip. 

The maximum bulk field was found to be independent of the beam angle 
whereas the peak tangential field at the edge naarest the beam source 
varied as where B is in degrees. The location of the maximum 

bulk field corresponds to the point where the electrons impact the surface 
at 0* and was found to shift coward the near edge for increasing beam 
angles. The charge distribution for non-normal incidence was also found 
to shift toward the edge nearest the electron source. For angles greater 
than 13'’ the maximum charge density was found at the neae edge. 

The surface potential and field profiles were found to be functions of 
Che normalized coordinate z/W over most of the sample surface , These 
quantities therefore can be described by a single scale length over a broad 
range of charging conditions. 

The external field quwtltles except at the edge were found to be 
largely independent of the sample thickness d. The internal bulk normal 
field was found to be inversely proportional to d. 

The simulation describing the charging of a specimen split by a central 
slot resembled that found for a uniform specimen over most of Its surface. 

A region of charge accumulation near the slot edges is strongly influenced 
by the relative slot width. For large slots the profiles should evolve 
toward those descriiiing two decoupled regions. External surface fields in 
the neighbourhood of the slot were found to be largely independent of the 
slot size and resembled the corresponding edge fields. 

The precharging simulations suggest that highly localized perturba- 
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tlons In the charge distribution could persist after most of the 
surface has attained a steady state. 

The nonuniform material simulations suggest large surface fields can 
develop in the neighbourhood of the thickness transition. 

The Kapton model, under emission-limited or conduction-limited 
charging conditions, produced two very different cross-sectional profiles. 

The conduction-limited cases produced a flat, uniform charge distribution 
whereas the emission-limited cases produced a broad central maximum. 

Maximum external surface fields were comparable although the tangential 
component associated with conduction-limited charging was very strongly 
localized at the edge. 

The Inclusion of tangential currents for Kapton had very little effect 
on the equilibrium state. Mo attempt was made to model a separate inherent 
surface conductivity, 

The internal front face field was found to be nearly constant for 
most simulations. Its magnitude was slgniflgantly smaller than the peak 
surface or bulk fields for both Teflon and Kapton under typical charging 


conditions. 
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